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SECTION  1 
INTRODUCTION 


1-1  BACKGROUND 

Systems  generated  electromagnetic  pulse  (SGEMP)  is  a  transient  radiation 
effect  produced  in  the  exposure  of  a  hardware  system  to  the  pulsed  x-radiation  and 
tf-radiation  from  a  nuclear  burst.  Electrons  are  driven  from  irradiated  material 
surfaces  which  generate  electromagnetic  fields  within  and  about  system  structure 
and  drive  currents  on  system  surfaces.  This  electromagnetic  environment  couples 
energy  into  electrical  systems;  if  sufficiently  intense,  this  coupled  energy  can 
produce  upset  or  damage  in  sensitive  electronic  components. 

For  satellite  systems  the  SGEMP  threat  can  be  classified  into  four  re¬ 
sponse  regimes:  external  SGEMP,  internal  SGEMP,  photon  direct  drive  SGEMP,  and 
black  box  SGEMP.  These  regimes  are  defined  in  the  following  paragraphs. 

External  SGEMP.  Electron  emission  from  irradiated  external  surfaces  of 
the  satellite  produces  a  space  current  which  drives  electric  and  magnetic  fields 
in  the  neighborhood  of  the  satellite  and  currents  on  structural  surfaces  of  the 
satellite.  These  currents  and  fields  may  couple  directly  to  electrical  systems 
through  coupling  into  antennas  and  exposed  electrical  cabling,  or  indirectly  through 
penetration  into  the  satellite  and  along  electrical  cabling  shields. 

Internal  SGEMP.  Photons  which  penetrate  into  the  satellite  interior  drive 
electrons  from  irradiated  internal  surface  producing  space  current  within  enclosures. 
These  cavity  space  currents  drive  electromagnetic  fields  and  replacement  currents 
which  couple  into  electrical  cabling. 

Photon  Direct  Drive  SGEMP.  Photons  which  penetrate  into  electrically 
shielded  cables  drive  electron  emission  from  the  interior  surfaces  of  cable  shields 
and  external  surfaces  of  cable  conductors.  This  charge  displacement  drives  replace¬ 
ment  currents  on  the  cable  conductors  which  propogate  to  interfacing  circuitry. 
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Black  Box  SGEMP.  Photons  which  penetrate  into  electrical  component  boxes 
drive  electron  emission  from  internal  surfaces  of  the  boxes  and  external  surfaces 
of  conductors  and  bulk  currents  within  dielectrics.  The  currents  drive  replacement 
currents  and  fields  which  couple  into  wiring  and  printed  circuitry.  This  effect  is 
particularly  significant  in  component  boxes  such  as  RF  amplifiers  which  utilize 
high  atomic  number  plating  materials  in  their  internal  construction. 

The  net  effect  of  the  various  SGEMP  responses  is  to  establish  an  electrical 
threat  to  the  interface  circuits  of  component  boxes  which  can  be  analytically  repre¬ 
sented  as  shown  in  Figure  1  by  a  number  of  separate  invariant  (independent  of  cable 
and  load  impedances)  distributed  drive  sources  The  field  coupled  threat  due  to  the 
external  and  internal  SGEMP  can  be  represented  by  capacitive  and  inductive  voltage 
sources,  Ec  and  EL  respectively.  In  the  absence  of  space  charge  limiting,  the  wave¬ 
shape  of  Ec  follows  the  photon  threat  pulse  shape  If,  while  E|_  follows  the  time  de¬ 
rivative  of  the  photon  pulse.  The  photon  direct  drive  threat,  on  the  other  hand, 
is  represented  by  a  current  source  whose  waveshape  al so  fol lows  the  photon  pulse.  It 
is  the  distributed  nature  of  the  threat,  the  voltage  or  current  source  representation 
and  the  specific  placement  in  the  transmission  line  circuit  which  gives  rise  to  the 
unique  energy-time  transmission  line  waveshape  profile  associated  with  each  SGEMP 
response  regime,  and  the  principal  reason  for  this  development  program. 


Figure  1.  General  SGEMP  wire-to-shield  drive  source  threats 
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It  is  possible  through  implementation  of  hardening  techniques  in  the 
satellite  design  to  strongly  reduce  the  susceptibility  of  the  system  to  the  exter¬ 
nal  SGEMP  to  a  level  small  compared  with  the  response  to  the  internal  SGEMP.  Ef¬ 
fective  techniques  are  (1)  the  containment  of  all  mission  critical  components  in 
Faraday  cage-design  enclosures  and  (2)  the  routing  of  all  cable  runs  between  enclo¬ 
sures  through  conduits. 

Similarly  the  system  susceptibility  to  internal  SGEMP  can  be  strongly 
reduced  through  implementation  of  hardened  design  techniques  to  a  level  small  com¬ 
pared  with  the  photon  direct  drive  SGEMP.  Effective  techniques  here  are  (1)  the 
use  of  RF-tight  design  for  all  mission  critical  component  boxes,  (2)  the  RF  shielding 
of  all  cables  and  connectors  interconnecting  mission  critical  components,  and  (3)  use 
of  low  Z  coating  on  all  internal  surfaces.  Therefore,  with  the  excention  of  black 
box  SGEMP  where  the  above  techniques  cannot  be  applied,  the  primary  concerns  of  the 
SGEMP  survivability  assessment  are  photon  direct  drive  SGEMP. 

1-2  PHOTON  DIRECT  DRIVE  THREATS 

The  consequence  of  the  photon  direct  drive  being  the  primary  SGEMP  threat 
is  that  the  threat  presented  to  interface  circuits  shown  in  Figure  1  is  altered  such 
that  the  capacitive  and  inductive  voltage  sources  become  negligible.  The  basic 
invariant  threat  becomes  a  distributed  current  source  (whose  waveshape  follows  the 
photon  pulse)  driving  a  wire  with  transmission  line  characteristics. 

Since  satellite  cables  are  comprised  of  many  unshielded  wires  enclosed  by 
a  common  bundle  shield,  the  actual  invariant  threat  presented  to  component  boxes  is 
as  shown  in  Figure  2  for  a  general  "N"  wire  cable.  This  threat  consists  of  "N"  in¬ 
dividual  distributed  current  sources;  one  for  each  wire,  driving  wi re-to-shield. 

The  physics  of  the  photon  direct  drive  interactions,  though,  are  such  that  large 
differential  mode  current  distributions  are  induced  in  the  cable,  that  is,  each  wire 
has  a  different  value  current  driver  induced  in  it.  This  mechanism  is  such  that  the 
outer  layer  wires  of  the  cable  core  on  the  X-ray  incident  side  have  high  current 
levels  induced  in  them.  The  source  currents  induced  in  the  wires  as  the  X-ray  pulse 
travels  through  the  core  become  progressi vely  smaller,  and  even  change  polarity 
sign.  Each  wire  in  the  cable  represents  a  real  wire-to-shield  transmission  line 
situation.  The  wires  on  the  outer  periphery  of  the  core  have  the  lowest  character¬ 
istic  impedance,  while  those  on  the  inner  layers  have  the  highest.  Each  of  the  "N" 
wires  on  the  bundles,  however,  are  also  tightly  coupled  to  one  another.  Hence,  the 
wires  in  the  cable  also  represent  a  real  wire-to-wire  transmission  line  situation. 
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Photon  direct  drive  threat  in  multiwire  satellite  cables 


In  spite  of  the  apparent  complexity  of  the  electrical  nature  of  the 
photon  direct  drive  effects,  the  threat  presented  to  each  interface  pin  can  be 
conceptually  reduced  to  an  equivalent  wire-to-shield  problem  and  represented  as 
shown  in  Figure  3.  The  source  term  current,  Iq^,  shown  here  is  an  equivalent 
source  term  which  represents  the  combination  of  the  actual  invariant  source  term 
induced  in  the  wire  and  the  amount  of  current  coupled  to  or  decoupled  from  the 
wire  of  interest  to  the  other  wires  in  the  bundle.  A  quantatative  representation 
sjch  as  this  for  each  wire  in  the  bundle  reduces  a  very  complex  problem  to  a  very 
simple  task  to  identify  which  wire  in  the  bundle  represents  the  maximum  energy 
threat  to  be  used  in  a  hardness  verification  test. 

The  above  discussions  have  alluded  to  the  importance  of  considering  the 
transmission  line  nature  of  the  photon  direct  drive  threat.  The  significance  of 
these  characteristics ,  which  will  be  discussed  in  full  in  Section  2,  are  implied 
in  the  typical  behavior  shown  in  Figure  4  and  5.  This  data  shows  the  magnitude  of 
the  threat  presented  to  component  box  interface  pins  in  terms  of  the  energy  level 
and  delivery  time  associated  with  a  1  ampere  per  meter  direct  drive  threat  in  a 
single  wire-to-shield  system  from  a  particular  photon  pulse  shape. 

It  will  be  shown  that  as  the  interface  impedance  value  increases,  the 
threat  energy  delivered  increases  until,  at  large  impedance  values,  the  threat 
capability  reaches  a  finite  maximum.  This  maximum  energy  capability  of  the  line 
increases  linearly  with  line  length.  However,  the  impedance  value  at  which  the 
threat  capability  saturates  decreases  as  line  length  increases  until,  the  threat 
is  essentially  a  constant  energy  value  independent  of  impedance  and  is  quite  simply 
Q^/2C  (i.e.:  the  equivalent  charge,  Q,  produced  by  the  X-ray  threat  divided  by 
twice  the  cable  capacitance,  C). 

It  will  also  be  shown  in  Section  2  that  the  time  it  takes  to  deliver  the 
threat  energy  is  also  markedly  effected  by  interface  impedance  .  For  an  interface 

impedance  which  is  matched  to  the  wire-to-shield  characteristic  impedance,  the  threat 
enemy  is  delivered  in  minimum  times  which  are  linearly  related  to  line  lenath.  Hiciher 

and  lower  impedance  values  result  in  longer  times  for  the  threat  eneroy  delivery. 

These  response  characteristics  demonstrate  the  complex  nature  of  the 
electrical  pulse  threat  delivered  to  interface  electronics  from  the  photon  direct 
drive  threat.  One  of  the  most  striking  features  of  the  response  is  that  the  load 
impedance  value  controls  the  amount  of  energy  it  receives  for  short  line  lengths. 


WIRE  NO.  N 


SHIELD 


COMPONENT  BOX  NO.2 - 1 


- COMPONENT  BOX  NO.  I 


Figure  3.  Equivalent  single  wire-to-shield  photon  direct 
drive  threat  for  component  box  interface  pins. 
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Figure  4.  Transmission  line  dependence  of  the  SGEMP  energy  delivered 
to  component  box  interfaces  from  a  1  ampere  per  meter 
photon  direct  drive  threat.  The  energy  values  shown  are 
50%  of  the  final  energy  delivered. 
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Figure  5.  Transmission  line  dependence  of  the  SGEMP  energy  rate 
delivered  to  component  box  interfaces  from  a  1  ampere 
per  meter  photon  direct  drive  threat.  The  time  value 
shown  is  the  time  in  which  50%  of  the  final  energy  is 
delivered. 


while  for  the  longer  lengths  the  threat  is  essentially  a  constant  energy  one.  The 
energy-time  characteristics  and  polarity  of  a  pulsed  electrical  threat  are  of  par¬ 
ticular  importance  in  determining  piece  part  burnout  vulnerability.  Hence  a  correct 
simulation  test  must  duplicate  the  actual  energy-time  nature  of  threat.  The  purpose 
of  this  work  was  to  develop  a  technique  and  instrumentation  to  drive  an  interface 
circuit  with  the  actual  waveshape  of  the  invarient  wi re-to-shiel d  source  term  which 
follows  the  photon  pulse,  and  with  the  actual  cable  length  and  characteristic 
impedance  of  the  cable  going  to  the  interface  circuit. 


1-3  MULTIWIRE  CABLE  CASE 

The  threat  energy  characteristics  given  in  Figure  4  show  that  the  maximum 
energy  capability  of  the  photon  direct  drive  is  limited,  as  would  be  expected,  to 
finite  values.  It  is  actually  the  wire-to-shield  characteristic  impedance  (which 
is  related  to  capacitance),  dielectric  constant,  and  wire  length  which  control  the 
SGEMP  threat  from  reaching  physically  unreal izeable  values.  The  maximum  energy 
capability,  in  fact,  is  inversely  proportional  to  dielectric  constant  and  directly 
proportional  to  characteristic  impedance,  line  length  and  total  charge  associated 
with  the  direct  drive  current  source.  This  surfaces  a  significant,  but  somewhat 
subtle,  aspect  of  the  photon  direct  drive  threat. 
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Recall  that  the  photon  direct  drive  threat  is  essentially  of  a  wire-to- 
shield  transmission  line  nature  with  large  differential  mode  source  terms  created 
in  the  cable  bundle.  Since  the  maximum  energy  capability  of  a  particular  wire  is 
not  only  controlled  by  the  magnitude  of  the  wire's  current  source  term  but  is  also 
controlled  by  the  wire-to-shield  characteristic  impedance,  it  does  not  necessarily 
follow  that  the  wire  with  the  largest  source  term  represents  the  maximum  threat 
condition,  since  this  wire  generally  has  the  lowest  wire-to-shield  characteristic 
impedance . 


In  reality,  each  wire  in  the  cable  has  it's  own  unique  maximum  energy 
capability  associated  with  it.  Conceptually,  one  could  describe  the  variation 
in  energy  threat  of  the  individual  wires  in  the  cable  core  as  a  function  of  wire 
characteristic  impedance,  as  shown  in  Figure  6  along  with  the  complete  mathematical 
expressions  for  maximum  wire-to-shield  energy,  short  circuit  current  and  open  circuit 
voltage. 


The  development  of  such  a  display  is,  in  fact,  what  would  be  required  in 
order  to  define  the  maximum  energy  threat  that  the  SGEMP  direct  drive  poses  to 
interface  circuits.  Furthermore,  the  definition  of  the  wire-to-shield  characteristic 
impedance  associated  with  the  maximum  energy  is  essential  to  correctly  simulate  the 
threat,  since  the  indiscriminate  use  of  a  lower  or  higher  impedance  value  would 
result  in  an  undertest  or  overtest  respectively. 

The  correct  threat  specification  situation  is  further  complicated  by  the 
fact  that  since  satellite  cables  are  comprised  of  many  unshielded  wires  enclosed  by 
a  common  bundle  shield,  significant  interwire  coupling  occurs  which  alters  the  equiva¬ 
lent  electrical  parameters  (equivalent  drive  source  level  and  wire-to-shield  character¬ 
istic  impedance)  of  the  transmission  line  type  threat.  From  a  threat  specification 
standpoint,  one  could  adopt  an  approach  based  on  defining  the  maximum  threat  posed 
to  component  box  pins  in  terms  of  the  total  cable  core  current  and  some  nominal  or 
worst  case  characteristic  impedance.  This  is  an  extremely  conservative  approach 
which  is  unrealistic  but  mathematically  representative  of  an  upper  limit.  This 
approach,  however,  is  completely  intolerable  from  a  satellite  system  standpoint 
since  it  represents  a  significant  hardware  impact  design  solution  to  providing 
adequate  hardening  against  a  realistic  threat  level. 
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Hence,  to  correctly  specify  and  simulate  the  maximum  threat  induced  by 
the  photon  direct  drive  response.  Sections  3  and  4  develop  analytical  techniques 
which  would  allow  one  to  define  a  realistic  set  of  electrical  parameters  which  are 
representative  of  the  realistic  maximum  energy  threat  contained  in  a  multi wire 
satellite  cable. 

1-4  REVIEW  OF  PRESENT  SGEMP  DIRECT  DRIVE  SIMULATION  METHODS 

Since  operational  testing  of  the  survivability  of  a  satellite  is  precluded 
by  the  nuclear  test  ban,  simulation  experiments  must  be  employed  to  perform  SGEMP 
hardness  verification  tests.  Photon  simulators  with  threat  level  fluences  are  not 
available,  however,  and  recourse  must  be  made  to  current  injection  test  (C.I.T.) 
techniques  are  required  by  the  present  satellite  system  survivability  specifications. 
The  various  C.I.T.  techniques  which  are  presently  available  under  existing  techno¬ 
logy  for  possible  use  in  the  hardness  verification  tests  are  as  follows: 

1.  Direct  pin  injection  from  a  commercial  pulse  generator 

2.  Inductive/Capacitive  (L/C)  coupling  to  cable  shield 

3.  Direct  injection  on  cable  shield 

An  advanced  development  technique,  which  ts  described  in  Section  5  of  this  report 
could  also  be  considered  as  "direct  injection  wire-to-shield" . 

The  primary  considerations  in  selecting  any  candidate  C.I.T.  technique 
are  associated  with  the  capability  to  reproduce  the  energy-time  characteristics 
of  the  photon  direct  drive  threat,  the  compatability  with  actual  satellite  equip¬ 
ments,  and  the  degree  of  risk  inherent  with  the  approach.  Detailed  evaluations  of 
each  candidate  approach  have  been  performed  by  GE/SD  to  identify  the  optimum  method 
which  should  be  used  to  perform  photon  direct  drive  simulations.  The  essential 
results  of  these  evaluations  are  shown  in  Table  1.  Direct  pin  injection  from  a 
single,  commercially  available  pulse  generator  is  not  included  in  the  Table.  This 
technique  could  conceivably  be  configured  to  provide  positive  current  control  at 
the  high  SGEMP  required  threat  frequencies  even  for  the  high  density  type  Ilarge 
number  of  connector  pins)  connectors  used  in  satellite  systems.  This  technique, 
however,  does  not  come  anywhere  near  being  capable  of  simulating  the  transmission 
line  nature  of  the  actual  SGEMP  threat.  The  conclusions  from  the  C.I.T.  evaluations 
are  that  the  C.I.T.  technology  which  has  been  developed  to  date  does  not  provide  the 
means  to  meet  the  SGEMP  hardness  verification  specification  in  design  hardened 
satellite  systems  where  the  SGEMP  threat  is  due  to  photon  direct  drive  effects. 


18 


Table  1.  Evaluation  of  various  current  injection  techniques 
for  photon  direct  drive  simulation 
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Techniques  such  as  those  previously  developed  under  EMP  programs  do  not 
provide  a  correct  simulation  of  the  energy-time  nature  of  the  photon  direct  drive 
threat  nor  do  they  provide  a  practical,  low  risk  approach  which  is  compatable  with 
photon  direct  drive  hardened  designs  in  satellite  equipments.  Capacitive  and  induc¬ 
tive  coupling  techniques,  for  example,  do  not  correctly  simulate  the  actual  trans- 
misson  line  nature  of  the  photon  direct  drive  source  term.  Furthermore,  these  tech¬ 
niques  do  not  provide  individual,  positive  control  of  currents  delivered  to  component 
box  pins  when  both  long  and  short  wire  threats  are  required  to  be  simultaneously 
injected  at  the  same  connector  interface. 

The  most  desirable  approach  identified  is  direct  injection  wire-to-shield. 
This  approach  is  the  most  desirable  from  the  standpoint  of  controlling  amplitudes 
on  individual  wires,  eliminating  unintentional  overstress  situations,  and  producing 
the  desired  waveshapes  to  simulate  the  SGEMP  direct  drive  threat.  The  technology 
development  effort  required  to  establish  a  realistic  C.I.T.  technique  which  is 
capable  of  both  correctly  simulating  the  electrical  nature  of  the  SGEMP  photon  direct 
drive  threat  with  a  multipin,  simultaneous  drive  capability,  as  well  as  meeting 
the  practical  and  hardened  desiqn  constraints  posed  by  actual  satellite  electronic 
and  complex  cable  harness  hardware  is  described  in  detail  in  Section  5. 

1-5  THREAT  SIMULATION  HARDWARE  CONSIDERATIONS 

Since  the  severity  of  the  SGEMP  photon  direct  drive  threat  is  controlled 
by  energy-time  transmission  line  characteristics  and  specific  wire  length,  hardening 
requirements  are  imposed  based  on  actual  cable  length  in  the  satellite  system  for 
each  interface  circuit  as  well  as  the  quantitatively  defineable  electrical  loads, 
not  only  at  the  circuit  of  concern,  but  at  the  interfacing  circuit  at  the  opposite 
end  of  the  wire.  In  many  cases,  for  example,  the  presence  of  a  particular  low 
impedance  electrical  load  or  the  addition  of  a  single  zener  diode  at  one  end  of  a 
short-to-medium  length  wire  is  sufficient  to  provide  protection  to  the  interfacing 
circuits  at  both  ends  of  the  line.  For  longer  length  lines,  SGEMP  hardening  is 
implemented  at  both  ends  of  the  line  consistent  with  the  line  length  of  the  trans¬ 
mission  line  type  threat.  Hence,  many  interfacing  circuits  in  a  component  box  are 
only  SGEMP  hardened  when  all  their  respective  interface  component  boxes  are  electri¬ 
cally  connnected. 
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An  additional  design  consideration  for  threat  simulation  is  associated 
with  the  physical  nature  of  the  electrical  systems  (cables  and  connectors)  in  actual 
satellites.  Satellite  cables  are  unlike  re-entry  vehicle  cables  in  that  they  are 
not  generally  straight,  two  connector  cables.  Instead,  satellite  cables  are  typical¬ 
ly  complex  configurations  comprised  of  large  bundles  of  upwards  of  hundreds  of 
unshielded  wires  contained  within  a  common  bundle  shield  with  numerous  branches  and 
connector  breakouts,  with  high  density  (large  number  of  pins)  connectors  ■  sed  for 
component  box  interfaces.  Furthermore,  the  wiring  in  these  cables  consist  of  wires 
with  short  total  lengths  (low  threat  level  wires  which  require  minimal  SGEMP  harden¬ 
ing  of  their  respective  interface  circuits)  in  a  common  bundle  with  wires  of  long 
total  lengths  (high  threat  level  wires  which  require  extensive  interface  circuit 
hardening) . 

The  consequence  of  these  design  aspects  are  that  an  SGEMP  simulation  tech¬ 
nique  must  not  only  provide  high  simulation  fidelity  of  the  energy-time  transmission 
line  nature  of  the  SGEMP  threat  for  multipin,  simultaneous  drive  of  component  boxes, 
but  must  also  be  capable  of  correctly  driving  a  component  box  whose  SGEMP  protection 
is  contained  in  numerous  other  component  boxes  which  are  not  in  the  test  setup,  as 
well  as  being  capable  of  simultaneously  driving  some  pins  at  high  levels  and  others 
at  low  levels  to  accomodate  various  distributions  of  long  and  short  wire  thrc'S  at 
the  same  connector  interface.  This  latter  requirement  also  has  a  significant  impact 
on  the  risk  potential  associated  with  any  potential  simulatin  '  techn^o  .*•. 

Since  the  SGEMP  hardness  verification  test  is  typically  performed  inline  at 
the  latter  portion  of  the  total  Qualification  Cycle  (generally  before  thermal -vacuum 
qualification),  a  low  risk  simulation  approach  is  required  to  avoid  inadvertent  over¬ 
stress  levels  which  could  cause  spacecraft  equipment  failures  and  impact  system  inte¬ 
gration  and,  possibly,  flight  schedules  In  view  of  the  typical  SGEMP  hardening 
designs  employed  in  satellite  systems,  one  cannot  generally  affort  the  luxury  of 
overdriving  some  pins  at  higher  levels  in  order  to  assure  that  the  minimum  required 
test  levels  are  properly  established  in  all  pins.  As  such,  the  simulation  method 
must  also  be  capable  of  providing  individual  positive  control  of  each  of  the  currents 
delivered  to  each  of  the  interface  pins  to  eliminate  high  current  (or  voltage)  hazards 
to  both  the  equipment  as  well  as  the  operating  personnel.  An  obvious  additional 
requirement,  which  is  more  practically  oriented  for  user  test  evaluation,  is  that  the 
simulation  methods  be  capable  of  being  configured  into  portable,  in-house  hardware 
readily  adaptable  to  correctly  interfacing  (both  electrically  and  mechanically)  with 
actual  satellite  component  boxes  and  their  associated  ground  test  equipments. 


SECTION  2 

TWO  WIRE  TRANSMISSION  LINE  CHARACTERISTICS 


2-1  BACKGROUND  AND  APPROACH 

This  section  describes  a  terminated  two  wire  transmission  line  which  has 
point  current  sources  distributed  along  its  length  from  return  to  signal  line.  A 
number  of  fundamental  transmission  line  response  equations  are  developed  which  pro¬ 
vide  useful  tools  for  analyzing  photon  direct  drive  cable  resopnses  in  a  X-ray 
radiation  environment. 

The  fundamental  equations  on  which  this  work  rests  appear  in  Schelkunoff 
(Reference  1)  where  the  waves  produced  on  a  transmission  line  by  point  sources  are 
succinctly  described.  The  work  of  King  (Reference  2)  and  Weeks  (Reference  3)  pro¬ 
vide  additional  development  of  transmission  line  theory  and  point  sources  which  form 
the  background  theory  to  this  section. 

The  response  of  a  loaded  two  wire  transmission  line  to  distributed  point 
sources  is  derived  in  Section  2-2.  The  point  current  sources  along  the  length  of 
line  are  assumed  to  be  proporational  to  a  defined  waveshape,  but  with  peak  amplitude 
and  start  time  arbitrarily  varying  along  the  length  of  line.  The  line  is  assumed  to 
be  linear,  time  invariant  (Reference  4)  and  lossy.  Loads  are  also  assumed  linear  and 
time  invariant.  The  equations  for  load  currents  and  for  current  on  the  line  are 
developed  with  reasonable  generality  in  the  frequency  domain. 

In  Section  2-3,  the  transmission  line  response  with  current  drivers  is 
then  studied  for  the  corollary  case  of  a  cable  illuminated  by  an  X-ray  with  an 
arbitrary  angle  of  incidence  so  that  the  induced  point  sources  turn  on  as  the  X-ray 
moves  along  the  cable.  Section  2-4  applies  the  results  of  Sections  2-2  and  2-3  to 
specific  applications  of  load  impedance,  pulse  width  and  a  current  source  amplitude. 

2-2  TRANSMISSION  LINE  RESPONSE  TO  DISTRIBUTED  CURRENT  SOURCES 

A  two  wire  transmission  line  of  length  S  is  shown  in  Figure  7.  The  line 
is  loaded  by  two  loads:  Z0  at  z=o  and  Zs  at  z=s,  where  z  is  the  distance  on  the 
line  from  its  end.  Both  line  and  loads  are  linear  and  time-invariant.  The 
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transmission  line  is  driven  by  a  distributed  current  source,  assumed  to  be  dis¬ 
tributed  point  sources  in  the  sense  of  Schelkunoff  (reference  1).  The  point  current 
source  is  described  mathematically  as: 

i0  ( x , t )  =  A(x)  g  (t-tx)  (2.1) 

where  A(x)  is  an  arbitary  amplitude  density  function  (Amp/m)  of  source  location  x, 
measured  from  the  end  of  the  line.  The  function  g(t)  is  assumed  to  be  dimensionless, 
normalized  and  causal,  i.e., 

[ g ( t ) jmax  =  1  (2.2a) 

and 

g(t)  =  0  for  t  <  0  (2.2b) 

The  parameter  tx  is  an  arbitary  function  of  source  location  x.  Thus,  the  equation 
for  i0(x,t)  is  a  representation  of  a  point  source  which  exists  along  the  length  of 
the  line  with  an  arbitary  amplitude  A  (x)  and  with  an  arbitrary  start  time  tx, 
dependent  on  location  x  from  the  end  of  the  line. 

The  Fourier  transform  of  i0(x,t)  can  now  be  written  as 


I0(x,w)  =  F  (i  (x  ,t) )  =  o(x^)  G(u) 

(2.3) 

where 

G(u>)  =  F  (g(t) )  =  f  g(t)  exp  (-jwt)dt 

and 

Jo 

p  (  X  ,  a)  )  =  A(x)  exp  (-ju)tx)  . 

For  simplicity. 

the  abbreviated  notation  I0(x,w)  =  I0(x)  and  p(x) 

is  introduced. 

2F 


Point  Current  Source  Density 
i0(x ,t)=A(x)g(t-tx) ,  Distributed 
Along  Line 

I(z)  =  Line  Current  Toward  Is  @  Location  I 
V(z)  =  Line  Voltage  @  Location  1 

Figure  7.  Transmission  line  with  distributed 
point  current  sources 


The  standard  transmission  line  notations  of  reference  2  are  introduced: 
characteristic  impedance  Zc  =  (2.4a) 

line  propogation  constant  =  V( R  +  ju>L)(G  +  jwC)  (2.4b) 

reflection  coefficients  ~  p  _  ^s  ~  ^c  (2.4c) 

zo  +  ?c  s  ls  +  2C 

where  R  =  series  line  resistance  P/M 
L  =  series  line  inductance  H/M 

G  =  shunt  line  conductance  S/M 

C  =  shunt  line  capacitance  F/M 


The  extensions  of  Zc  and  ^  to  include  skin  effects  and  other  high  frequency  losses 
are  covered  in  reference  2. 

The  objective  of  this  derivation  is  to  compute  the  current  on  the  line  Iz 
at  location  z.  The  current  Iz  can  be  thought  of  as  the  difference  of  two  currents: 

I*  =  I l (2)  -  I2(Z)  (2.5) 

where  I](Z)  is  the  contribution  to  I2  of  all  sources  to  the  left  of  z, 
i .e  ,  x<z  and  1 2 ( Z )  is  the  contribution  of  all  sources  to  the  right  of  z,  i.e., 

X>2  . 


Using  Schel kunoff ' s  development  (reference  1),  a  point  source  driver  at 
location  x  of  Zs  to  location  x  and  where  Z|_  is  the  reflection  of  Z0  to  location  x 
as  indicated  in  Figure  8A.  It  is  assumed  that  a  differential  source  term  driver 
I0(x)  dx  within  a  differential  length  of  line  dx  launches  a  wave  dlR(x)  with  the 
sign  convention  of  Figure  8a  indicating  plus  to  the  right.  Thus,  by  current  division. 


where 


Zl>) 

dlR(x)  =  ZL(x)  +  Zr(x) 


1  +  'o  exp  1 

;-2/x 

1 

1  -  p  exp  l 
1  +  Q  exp  1 

l-2j 

t-2> 

f  x 

t 

s-x)) 

1  -  p  exp  1 

(-2: 

FT 

5-x|) 

Reflecting  dl r ( x )  to  location  1  gives: 


dl- 


exp  1 

*)*  's 

A 

M 

exp  1 

RxT-  p  exp 

M 

^2?-x)T 

(2.6a) 

(2.6b) 

(2.6c) 


(2.7) 


Rr  ( ^ » x )  dIR 

where  differential  operator  "d"  is  on  x  only. 


Using  the  complementary  model  of  Figure  8b  for  the  left  current  (sources  at  x>2) 
yields 


zr(x) 

(2.8) 

dIL(x)  =  znTjnrziTxT  Io(x)  dx 

and 

(2.9) 

r  Rp2,x)dlL 

•>  7 


a)  Location  x<z. 
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?,(*)  y  **(<) 

I(*)d.x 


Figure  8a. 
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b)  Location  x>z.  i . 


— <  r 
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Figure  8b. 


Then  from  Equation  2.5  I j  is  written  as: 


A  ■  f 


The  integrals  are  of  the  form: 


|  Hr  (Z,X)  Iq(X)  dx  -  |  Hl  (Z,X)  Iq(X)  dx 


where 


zi  W 

hr  =  Rr(Z,X)  zl(x)+zr(X) 

ZR  (  X  ) 

HL  =  Rl(Z’X)  z.  (X)+Zn(X) 


Contrary  to  the  devastating  appearance  of  these  expressions,  term 
cancellations  occur  throughout,  so  that  Equation  2.12  is  all  that  remains. 
One  is  careful  to  note  that  the  direction  of  Iz  is  assumed  to  be  towards  Zs. 


e->(2w)] 

2( 1-r  r  e'2<s) 


[  e'2-r„  e-’2' 


eYX+r  e-,x 
0 


p(X)dx 


!(1-r» 


■  n 


e">x-t-rs  e‘2>s  e|X  j  p(X)dx 


(2.12) 


Equation  2.12  is  the  desired  transfer  function  Iz/G(w)  which  relates  Iz 
to  the  input  source  pulse  shape  G(w).  The  integrands  are  of  simple  exponential 
formats,  with  p(x)  =  A  (x)  exp  (-jwtx).  The  Fourier  spectrum  for  line  current  Iz 
is  thus  easily  computed,  once  p(x)  is  defined.  Time  domain  histories  can  be 
obtained  via  standard  Fourier  inversion  for  a  specified  G(w). 

The  line  voltage  at  location  z  can  now  be  computed  from  the  line  current 
Iz  via  superposition  from  the  relation: 


v ( z )  =  ZR(Z)  J  dij  -  ZL(Z)  | 


(2.13) 


where  the  functions  Zr  and  have  been  identified  in  Equations  2.6b  and  2.6c. 
The  integrals  for  I]  and  are  easily  identified  from  Equations  2.10  and  2.12. 

The  load  current  in  Z$,  namely  I^,  is  imply  Iz  evaluated  at  Z=S,  so 


IL  u 
L2 


e'rS(l-:'  ) 


^■Vs6"'5)  Jo 


rs 


e*,x  P(X)dx 


(2.14) 


The  load  current  into  Z0,  namely  I[_i  with  the  sign  convention  of  Figure  7  is  -  I. 
evaluated  at  z=0. 


'L1 U  0--n>  f 

°U)  2(‘-:o  's  e"2>$)  'o 


» 

^e",x+rs  e"2>s  e'x  j.(X)dx  (2.15) 
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2-3  APPLICATION  TO  X-RAY  DIRECT  DRIVE  RESPONSE 


« 


The  application  of  Equation  2.12  to  cable  direct  drive  resulting  from 
X-ray  illumination  is  straight  forward.  The  geometry  of  Figure  9  shows  a  cable 
with  an  X-ray  illumination  environment.  The  distributed  current  source  is  assumed 
to  have  the  form  of  equation  2.1: 

i0(x,t)  •  i„(t-tx>  =  I0  g(t-tx)  (2-16) 

where,  A(x)  =  ID  =  constant.  If  the  angle  of  incident  illumination  is  0,  as  shown 
in  Figure  9,  the  point  drivers  turn  on  along  the  line  at  time  tx  =  x  coso/vo  where 
V0  is  the  speed  of  light  in  M/S  and  tx  is  time  relative  to  the  time  point  when 
location  x=o  is  illuminated. 

By  taking  the  Fourier  transform  of  i0,  the  function  of  p(x)  is  identified 
as: 

p(x)  =  Ip  exp  ( -ju>X  cose/v0) 

(2.17) 

=  ID  exp  (-jk<0  cosox) 

Thus,  for  normal  incidence,  tx=o  and  the  entire  cable  is  driven  simultaneously, 
resulting  in  p(x)  =  IQ  =  constant.  For  on  line  incidence,  h  =  o  and  tx  =  x/vq, 
resulting  in  p(x)  =  IQ  exp  (-Jujx/vo). 


i 


Figure  9.  Cable  geometry  with  radiation/ 
incident  at  arbitrary  angle  0-  0-  ^ 
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In  other  words,  when  the  fluence  propagates  along  the  length  of  cable, 
the  X-ray  illuminates  the  cable  with  sources  that  turn  on  at  time  x/vQ.  The  cable 
generates  a  wave  which  moves,  generally  speaking,  slower  than  v0.  Equations  2.14 
and  2.15  will  now  be  studied  to  find  out  what  effects  these  wave  motions  have  on 
load  currents. 


Whereas,  the  equations  for  load  currents  I|j  and  1 1_2  (2.14  and  2.15)  are 
perfectly  general  and  can  be  analyzed  via  numerical  Fourier  analysis  for  lossy 
cables,  the  consideration  of  lossless  transmission  lines  is  now  taken  up.  This 
restriction  is  applicable  for  the  analysis  of  many  practical  cable  configurations. 
The  lossless  line  also  allows  for  the  introduction  of  simple  line  parameters,  in 
particular,  the  line's  one  way  transit  time  t,  so  as  to  provide  a  continuing 
physical  insight  to  the  meaning  of  the  equations  as  they  are  developed.  The  follow¬ 
ing  notation  is  introduced  (R=G=0): 


vo  = 

3  x  10^  M/S,  speed  of  light 

(2.18a) 

S  = 

length  of  line  in  meters 

(2.18b) 

To  = 

S/vQ  =  time  for  light  to  traverse  line  length  S. 

(2.18c) 

e0  = 

c/v0 

(2. 18d) 

vp  " 

( 2 . 1 8e ) 

T  = 

S/v  (line  delay  or  one-way  cable  transit  time 

(2. 18  f ) 

^  of  a  wave  on  the  line.) 

Y  = 

j^vTc  =  je 

(2. 18g) 

Ti  3 

H 

o 

o 

o 

</> 

(2.1 8h ) 

l2  = 

t+to  COS" 

(2. 18i ) 

3 

F(V*>) 

(2.1 8 j ) 

SI-C-) 

Q(w)  =■ 

0 

(2.18k) 

The  result  of  substituting  ,  (x)  =  exp  (-jr;  cos^  x)  into  the  Equation  2.14 


I|_2  yields  Equation  2.19,  which  when  factored  in  yields  Equation  3.20. 


for 

These 


equations  are  valid  for  i  >  t  cosh. 
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( l-r  )  Q(w)  <, 

!L2  =  2 [1  -rQi  s  exp  (-2j^T)]|  rj  exp  ("juJTo  cose) 


--f2.exp[-j,(2t>To  cos.)]  -(^-'T^exp 


(  -  jco  T  ) 


-  1  -r0rs  exP  (-2J-0 


(1  -•  s )  Q(u)  -J-  0  cos 
lL2  =  TTjT  e 


(2.19a) 


(2.19b) 


l'-4 


(2.20) 


Considerable  algebraic  effort  is  saved  in  the  integral  evaluation  for  IL] 
by  using  the  solution  for  I j_2  with  Zo  replaced  by  Zs  in  Figure  9  and  noting  that  in 


this  geometry. 


(X)  =  exp  (-ji-  (S-X)  cos. 


(2.21) 


Thus,  to  obtain  IL1  from  Equation  2.20  for  1^,  is  replaced  by  -t^,  Zo  is 
interchanged  with  Zslthereby  T  is  interchanged  with  -:  ) ,  j  is  interchanged  with 
l2’  To  is  rePlaced  bX  'T0  and  the  whoie  expression  is  then  multiplied  by 
exp  (-jii  S  cos'.)  =  exp  (-jwiQ  cos  •)  to  yield: 


(1-rJQU) 


•(■•  — K" 2  i  <*•< 


with  A  defined  in  Equation  2.19b.  The  student  of  Fourier  theory  will  immediately 
recognize  that,  for  o/o ,  will  start  responding  before  1^,  because  of  the 

absence  of  the  phase  shift  multiplier  exp  (-jo/'  cos  o)  which  occurs  in  Equation 
2.20.  That  is  to  say  that  the  current  in  the  load  at  z=s  will  be  zero  for  time 
less  than  (r  cos  o),  the  time  at  which  the  X-ray  arrives  at  z=s,  whereas,  the 
current  at  z=o  will  commence  upon  the  illumination  of  that  end  of  the  cable. 
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For  the  case  of  perpendicular  incidence,  where  the  line  is  instantaneously  and 
uniformly  illuminated,  o  =  -n/2  and  t=Tj=t2  yield: 


(1-r  )Q(u> ) 


1-r  e 
s 


-2j"T  -(l-rs)  e'ju)Tj 


(2.23) 


(1-rjQd) 


l-i0  e‘2ju)X  -(l-i 0)  e'ju)T 


where  A  =  1  -*'0i  s  e*P  (-2d  i)  as  defined  in  Equation  2.19b. 


(2.24) 


The  term  Qd)  has  been  defined  in  Equation  2.18k.  This  function  is 
identified  as  the  total  charge  accreted  on  the  open  circuited  line. 

Qd)  =  F  (q(t)j  (2.25) 


where 


q(t)  -  S  I  dd)<k 


(2.26) 


It  is  interesting  to  note  that  the  load  currents  are  proportional  to  line 
charge  accretion,  and  that  the  cable  in  a  sense  tries  to  integrate  the  waveshape  of 
the  illuminating  X-ray.  A  detailed  exposition  of  the  wave  structure  follows. 

2-3.1  Time  Domain  Response  of  a  Resistively  Loaded  Lossless  Line 

The  frequency  domain  equations  for  I^  and  1 ^  can  be  inverted  to  the 
time  domain  by  Fourier  theory.  The  equations  for  load  currents  are  of  the  form: 


l  -  Q  d)-= _ Q  1^)  ■  rp - * 

L  A  1  -"0rS  exp  ("2joJT' 


Using  the  binomial  series  expansion  for  1/.:  (reference  2),  1^  becomes: 

1 L  =  Q  d)  £  O'0r$)n  exp 


(2.27) 


(2.281 


so  that  for  resistive  loads,  i.e.,  r  and  r_  both  real, 

o  S 


1,  (t)  =  f  (r  rjn  f (t-2nx 


(2.29) 
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where  f(t)  is  causal  so  that  the  upper  limit  of  the  suirmation  of  Eauation  2.29 
may  be  changed  from  <*=>  to  n  ,  and  the  following  definition  is  implied: 


(r0rs)n  =  1  for  n=0 
when  either  rQ  or  r  or  both  are  zero. 

Noting  too,  that 

Q(w)  exp  (-jwt0)  =  F  (q(t-tQ)) 

enables  one  to  write,  for  the  lossless  resistively  terminated  line: 

i.  (t) 

L1 

I 


L2  "  2t1  (nto  lr0rS)  (t'To  cos9  -2nx) 

-  1  (rorS>"  (‘"<2"+1>') 

-  "f  (r0rs)"(-^)q  (t-x0  cose  -<2r*2|x)  | 
where,  t  has  been  defined  in  Equation  2.18: 

These  equations  are  valid  only  for  t>t0  fcoso).  The  case  of  an  air  line,  i .e 
is  addressed  separately  below. 

The  currents  in  either  resistive  load  are  sums  of  the  waveforms  as 
in  Equation  2.26  for  tiO  and  g(t)=0  for  t<o  for  the  accreted  charge. 


(1-r  )  { n-c° 

—  !  I  (r „rc)n  q(t-2m) 


2'2 


n=°o  /  1 

-  2  (rQrs)n  (l-rs— jq  (t-xo  cose  -(2n+l)T) 


n=o 

n=oo 

I 

n=o 


l  (rorS>"(-^->  (‘-<2-2> 


(l-rs)  (n.. 


(2.30) 


(2.31) 


(2.32) 


(2.33) 


■  >  f  =  T0 

given 
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These  summations  can  be  regrouped  by  introducing  two  functions,  which  are  dimension- 
ally  currents: 


V1)  =  ^  "  q^t-To  COS0"T) 

ib(t)  5  27J”|q^t-To  COSO) 


(2,34a) 

(2.34b) 


The  load  currents  become: 


n=° 


Vj")  ■  ‘Vs>"  'a 


n=o 


n=c 


*  <1-ro)rS  2  (rorS>  <b 
n=o 


n-00 


i.  (t)  =  (1-r, )  2  (r  r<-)n  i.  (t-2m) 
l2  b  n=o  0  b  b 


n-°° 


+  (1_rS)ro  ?  (rorS)n  ia  (t-T-2nT) 


n=o 


(2.35) 


(2.36) 


The  structure  of  Equation  2.35  allows  one  to  see  that  the  initially 
incident  pulse  on  load  Zo  is  ia.  This  pulse  is  followed  by  re-reflectance  of  ''a 
off  load  Zs,  plus  the  reflection  of  ’b  initially  incident  on  load  Zs,  plus  its 
re-reflections . 


2-3.2  Resistively  Matched  Line 

Equations  2.35  and  2.36  are  vastly  simplified  when  the  line  is  matched 
at  both  ends,  i.e.,  when  r0  =  ;s  =  0.  In  this  case,  there  is  no  reflection  off 
either  load,  and  the  load  currents  are  equal  to  the  first  incident  pulse: 

i^U)  =  i a ( t )  =  -2“  jq(t)  'q(t-T0  C0S°-T)j  (2.37) 

iL2(t)  =  ib(t)  =T^“|q(t'To  COSO)  "q(t-T)|  (2-38) 

again  subject  to  the  restriction  that  t  >  1  coso. 


2-3.3  The  Electrically  Short  Cable 


The  case  of  a  uniformly  illuminated  short  cable  provides  some  additional 
clarification  of  the  meaning  of  Equations  2.35  and  2.36  for  load  currents.  For 
simplicity,  the  case  where  0=11/2  is  considered,  for  which  (Equation  2.34)  the 
incident  pulses  become: 


where 


,  (t )  =  i h(t)  = 


-^-|q(t)  - q ( t-T  )| 


q(t)  =  S 


d^  with  A(x)  =  constant 


(2.39) 


From  Equation  2.35,  the  current  in  load  1  becomes: 


i,  (t) 
L1 


(1 


-ro>  2  <ro!'s>"  I  ia(t-2" 
n=o  l. 


t)  +  rs i a ( t-x -2riT  )  >  (2.40) 


The  current  in  load  2  is  also  given  by  Equation  2.40  with  r0  interchanged  with  rs, 
because  the  cable  is  assumed  to  be  instaneously  and  uniformly  illuminated. 


Since  the  one  way  cable  transit  time  is  defined  by  Equation  2.18f,  the 
pulse  ia  strictly  goes  to  zero  as  t  goes  to  zero,  and  therefore,  iL1  goes  to  zero 
too.  However,  if  one-  writes  for  i a ( t )  the  integral  expression 

ia(t)  =Triq(t)  '  q(t~T)  l  =  ir  J  (2-41) 

t-T 

one  may  approximate  the  integral  by 


|  io(^)d;  -  T  io(t) 


(2.42) 


if  i  is  slowly  varying  over  the  time  interval  t.  This  relation  provides  adefini- 

tion  of  what  constitutes  an  electrically  short  cable  with  respect  to  a  drive  pulse 

i0(t).  With  this  relation,  i  (t)  simply  becomes: 

a 

ia(t)  =  I  io(t)  (2.43) 
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Now  if  load  1  is  matched  so  that  rQ  =  0  and  load  2  is  open  so  that 
r  =  1  ,  then  from  Equation  2.40 


i|_  (t)  =  iaU)  +  ia(t-x) 


(2.44) 


independent  of  cable  length.  When  xis  small  so  that 


ia(t)  =  ia(t-x)  =  |  i0(t) 


(2.45) 


then  the  well  known  result 


iL1(t)  =  S  io(t) 


(2.46) 


is  obtained. 

The  case  where  load  1  and  load  2  are  equal  and  resistive  is  somewhat  more 
complex,  but  tractable.  In  this  case  rQ  =  r<.  =  (r|-l  )/(q+l  )  where  =  R/Zc .  Defining 
r  =  r0  -  l’s  allows  one  to  write  from  equation  2.40. 


i,  (t)  =  2 
LI 


i  (t-2nt)+r 

a 


( t-T-2ni ) 


(2.47) 


again  independent  of  cable  length.  Now  assuming,  for  any  value  of  n,  that  i  is 
small  enough  so  that 

i  (t-2m )  =  i  (t-i -2nx )  (2.48) 

a  a 


al lows  one  to  wri te 


1L1lt)aZ^Tnljr2n(1+r)ia<t'2n,)l 

One  may,  with  suitable  caution,  go  further  and  assume  that 
i  a ( t-2nx )  =  i a  ( t ) 


(2.49) 


(2.50) 
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for  all  times  of  interest,  and  then  by  identifying  the  binomial  expansion  for 
2-1 

(1-r  )  one  obtains: 

vt)'2W(-nr),*(t,-,.<t»  <2'5i) 

With  i  (t)  =  S/2  i  (t)  one  finally  obtains  for  the  short  cable 

iL  (t)  =  |  iQ(t)  (2.52) 

This  result  simply  states  that  the  cable  source  term  divides  equally 
between  the  two  loads.  However,  the  application  of  Equation  2.52  to  any  specific 
case  depends  on  load  resistance  R  and  pulse  shape  i  .  Clearly,  for  a  rapidly  vary¬ 
ing  i  (t),  or  a  bipolar  i q ( t ) ,  both  2.42  and  2.50  may  be  readily  violated,  and 
effects  arising  from  cable  energy  storage  may  become  important  to  the  computation 
of  load  current  rise  time/fall  time,  as  well  as  to  the  duration  of  the  trailing 

edge  of  the  load  current  pulse.  The  effects  of  r  are  also  implicit  in  the  use  of 
2N 

2.50.  As  r  approaches  an  acceptably  small  number  so  that  Equation  2.50  need  only 
be  valid  for  n-N,  then  2.52  is  an  acceptable  approximation.  Therefore,  the  closer 
r  is  to  zero,  the  less  terms  will  be  necessary  in  2.49  for  convergence  and  the  more 
accurate  2.52  will  be  for  a  given  i  (t)  and  given  t  small.  But  even  for  small  -, 
when  |r|-»l,  Equation  2.52  is  not  valid.  As  r  * +1  ,  the  cable  capacitance  determines 
the  load  pulse  shape,  since  the  cable  will  begin  to  charge  faster  than  charge  can 
flow  thru  R.  As  r  ► -1 ,  the  cable  will  ring,  albeit  at  a  high  freuqency,  because 
the  energy  absorbed  per  cycle  is  low,  regardless  of  x.  Clearly,  when  r  =  -1  ,  the 
cable  is  truly  shorted  and  being  lossless,  will  ring  forever  so  that  2.52  is  not 
valid  at  all. 

Thus,  the  derivation  of  Equation  2.50  is  offered  as  a  heuristic  indication 
that  the  more  complex  expression  of  Equation  2.40  is  consistent  with  elementary 
Kirchoff's  laws  applied  to  lumped  element  representations  of  the  cable  source 
problem. 
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2-4  ILLUSTRATE  GENERIC  RESULTS 

Because  the  load  current  pulse  waveforms  previously  derived  are  quite  com¬ 
plex  summations  of  pulses,  this  section  provides  some  numerical  and  graphical  re¬ 
sults  which  highlight  the  nature  of  cable  direct  drive. 

In  Section  2-4.1,  important  energy  relations  are  provided  which  relate 
the  absorbed  load  energy  to  load  resistance,  to  time,  and  to  cable  length. 

In  Section  2-4.2,  the  effects  of  angle  of  incidence  and  cable  length  for 
a  square  pulse  driver  on  a  matched  line  are  developed. 

2-4.1  Energy-Time  Relations  in  Various  Resistive  Loads 

An  important  phenomenon  in  a  cable  with  distributed  current  sources  is  that 
the  energy  which  can  be  extracted  from  the  cable  is  a  strong  function  of  the  termina¬ 
tion  loads.  This  section  provides  some  insight  as  to  what  these  terminal  effects  are. 

Equations  2.23  and  2.24  have  been  analyzed  via  a  numerical  Fourier  inverse 
with  a  uniformly  distributed  critically  damped  current  driver  of  the  form: 

1o(t)  =  Io(~to)exp  (>-4r)u(t>  (2-53) 

where  Ip  =  peak  current  source  in  amps  per  meter,  tm  =  time  of  peak  amplitude,  and 
u(t)  =  Heaviside  unit  step  function.  The  cable  is  assumed  to  be  instantaneously 
illuminated,  so  that  o  =  *72.  The  cable  is  assumed  to  be  loaded  at  either  end 
with  a  resistor  R  so  that  Z0  =  Zs  =  R  and  r0  =  rs  =  r.  With  Ip  normalized  to 
unity,  the  energy  in  each  load  has  been  obtained  for  study  by  a  computation  of  the 
running  integral 

EL(t)  ■  R  jVkld;  (2-54> 

0 

where  ip(t)  is  the  load  current  as  a  function  of  time.  The  integration  was  per¬ 
formed  via  a  simple  Simpson's  rule  from  the  time  trace  obtained  by  the  numerical 
Fourier  inversion  of  the  current  spectrum. 
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The  results  of  these  computations  are  shown  in  Figure  10  where  the 
normalized  energy  has  been  plotted  as  a  function  of  e  =  R/Zc  with  t  =  <*>.  E|_  has 
been  arbitrarily  normalized  to  indicate  relative  behavior  of  Eg  versus  5  and  cable 
length  S.  Since  S  is  proportional  to  cable  transit  time  t  (Equation  2.18f),  S  has 
been  indicated  by  i.  The  time  Tgg  is  a  measure  of  driver  pulse  duration,  and  is 
equal  to  the  time  when  the  driver  pulse  delivers  90%  of  its  total  energy  to  a  re¬ 
sistive  load  via  the  action  integral  relation: 

r  ?  rT9o  o 

0.9  I  iQ2(t)dt  =  I  iQ2(t)dt  (2.55) 

0  0 

The  results  of  Figure  10  show  that,  as  cable  length  is  increased,  the 
total  energy  absorbed  by  the  loads  is  increased.  Furthermore,  Figure  10  shows 
that  as  load  impedance  is  increased,  the  energy  the  load  absorbs  is  also  increased. 
For  the  geometry  analyzed,  equal  loads  at  either  end,  each  load  receives  one-half 
the  energy  delivered  by  the  cable.  In  the  limit,  for  R>>ZC,  the  loads  absorb  all 
the  energy  in  the  cable  when  open  circuited,  namely: 

E»ax  ■  I  %2/C  (2'561 

where  Qg  is  the  total  accreted  charge  and  C  =  t/Zc  is  the  total  cable  capacitance. 
For  the  pulse  shape  under  consideration  Qg  is  given  by: 

Q0  .  S.„  .  «.  ■  £  q(t) 

where  e  is  the  Napierian  base.  Since  Qg  is  proportional  to  cable  length,  as  is 
total  cable  capacitance,  the  maximum  energy  available  from  the  cable  is  also  pro¬ 
portional  to  cable  length. 

On  the  other  hand,  for  R<<ZC,  the  energy  absorbed  is  significantly  less 
than  the  total  energy  available  in  an  open  circuited  cable.  The  physical  explana¬ 
tion  of  this  phenomenon  lies  in  the  fact  that  low  impedance  loads  suppress  the 
line  voltage  as  incident  waves  are  reflected,  thereby  reducing  the  amount  of 
energy  the  cable  can  provide.  It  is  an  important,  as  well  as  interesting  observa¬ 
tion  to  note  that  the  current  drivers  are  assumed  independent  of  line  voltage  and 
that  consequently,  the  conservation  of  charge  flow  is  the  controlling  law,  making 
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zed  load  enerqy  versus  normalized  load  resistance 


the  available  cable  energy  strictly  a  function  of  the  voltage  level  which  the 
cable  is  allowed  to  build  up  to. 

Of  interest  in  Figure  10  is  the  knee  in  the  energy  curve,  i.e.,  the 
value  of  (R/Zc)  below  which  energy  falls  off  as  (R/Zc)  is  decreased.  The  location 
of  the  knee  is  a  function  of  the  electrical  cable  length,  relative  to  drive  pulse 
duration  Tgg.  A  cable  possessing  a  one  way  transit  time  t  which  is  signficantly 
longer  than  the  duration  of  the  drive  pulse  i0  will  almost  fully  charge  before 
the  effects  of  end  terminations  can  act  to  reduce  line  voltage.  Thus,  longer 
cables  are  capable  of  delivering  higher  energies  regardless  of  termination. 

Some  specific  examples  have  been  plotted  in  Figures  11  through  15.  The 
current  driver  has  been  taken  as  a  square  pulse  with  a  constant  Ig’Tp  product  of 
30  Ans/m.  For  these  figures  the  cable  length  has  been  held  constant  at  3m  with 
a  dielectric  constant  equal  to  9.  The  curves  were  derived  from  Equation  2.54 
using  Equations  2.32  and  2.33  as  the  load  currents  through  pure  resistive  loads. 
The  five  plots  present  a  family  of  curves  for  loads  between  0.0 5ft  to  50kn  and  for 
pulse  widths  between  300ns  to  3ns. 

The  effect  of  the  dielectric  constant  on  the  rate  o^  energy  absorbed  in 
the  load  resistance  for  values  of  Z|_  between  0.05si  to  50ks2  is  shown  in  Figures 
16  through  19  for  value  of  e=l  (vacuum)  to  81  (water).  For  values  of  most  typical 
cable  dielectric  materials  with  er  between  1  and  3  the  variation  in  energy  for 
matched  loads  is  relatively  small  and  increase  to  nearly  a  factor  of  two  between 
these  limits  for  values  of  r^>>!  or  <<1. 

The  effect  of  the  cable  length  on  the  temporal  dependence  of  the  energy 
absorbed  in  the  load  resistance  is  shown  in  Figures  20  through  23.  As  the  cable 
length  increases  from  0.3m  to  30m 

E ( t ,t|) — ►E(t,n”l ) 

W’>>-  Emax<ZL=Zc> 

That  is  the  maximum  energy  in  the  load  impedance  becomes  independent  of  the  load 
impedance  for  long  cables.  Furthermore,  the  maximum  energy  limit  increases  as  the 
cable  length  increases  from  0.033u.)  for  S  =  0.3m  to  33uJ  for  S  =  30m. 
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The  effect  of  the  characteristic  impedance  on  the  absorbed  energy  is 
shown  in  Figure  24  for  Zc  =  5ft  and  500ft  and  a  family  of  load  resistances  between 
0.005s:  and  500kft. 

The  energy  limit  (Emax)  which  is  absorbed  in  the  load  impedance  is  sum¬ 
marized  in  Figures  25  to  28  as  a  function  of  s,  the  load  impedance  normalized  to 
the  cable  characteristic  impedance.  The  summaries  show  the  variation  of  Emax 
versus  n  in  terms  of  a  set  of  parametric  curves.  The  parameter  variables  modeled 
are:  1)  pulse  width  (TD=3  to  300ns)  with  constant  total  charge,  i.e.,  I0'Tp  = 
constant  (30A'ns/m)  in  Figure  25,  2)  dielectric  constant  (e=1  to  81)  in  Figure  26, 

3)  cable  length  (s=0.3  to  30m)  for  three  values  of  er  (1,3,9)  in  Figure  27  and 

4)  cable  impedance  (Zc=5  to  500ft)  in  Figure  28. 

Figures  29  and  30  summarize  the  maximum  energy  absorbed  in  the  load  as 
a  function  of  the  cable  dielectric  constant  for  a  family  of  normalized  load  impe¬ 
dances.  The  two  figures  consider  the  necial  cases  for  cable  characteri Stic 
impedance  of  50ft  and  free  space  (377.:).  Figure  31  shows  the  variation  of  absorbed 
energy  with  dielectric  constant  similar  to  Figure  29  but  with  the  pulse  width 
increased  from  30ns  to  90ns  for  the  same  I0  (lA/m).  In  general  the  absorbed 
energy  increases  as  the  pulse  width  increases  at  a  rate  which  is  higher  for  larger 
n.  An  analogous  situation  is  illustrated  in  Figure  32  where  the  cable  length  is 
shortened  from  3m  to  0.3m  in  comparison  to  the  parameters  in  Figure  29  with  the 
pulse  width  held  at  30ns. 

Equally  important  as  the  total  energy  into  each  load  is  the  rate  at 
which  the  energy  is  deposited,  since  energy  rate  is  a  measure  of  dissipated  power. 
As  indicated  by  the  time  expansion  for  i |_ ( t )  with  n=  in  Equation  2.40,  the  load 
current  is  a  series  of  pulses  whose  shape  depends  on  line  delay  t  and  reflection 
coeffi cients  P.  Thus,  the  instantaneous  power  in  each  load  can  be  expected  to  peak 
periodically.  In  a  real  system,  the  instantaneous  power  may  not  be  as  important 
as  the  total  energy  accrued  in  the  load  by  N  successive  reflections.  In  order  to 
measure  the  energy  rate  then,  an  arbitrary  time  reference  is  introduced,  namely 
the  time  at  which  90  percent  of  the  total  energy  is  absorbed  by  that  load,  denoted 
tgg.  The  underlying  reason  for  this  time  reference  lies  in  the  fact  that,  theoreti¬ 
cally,  a  mismatched  cable  will  have  finite,  albiet  vanishingly  small  currents  for 


56 


ENERCY  -  /iJ 


Figure  26.  Maximum  absorbed  load  energy  vers 
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absorbed  load  energy  versus  normalized  load  impedance  for  various  line  impedances 


times  much  greater  than  t 


gQ.  Mathemeti cally ,  tgQ  is  obtained  from  the  relationship: 

EL(t9Q)  =  0.9ElH  (2.58a) 

or  via  the  equivalent  action  integral  relation 

90  9  r  o 

iLZ(t)dt  =  0.9  1  (t)dt  (2.58b) 

o 

A  plot  of  normalized  tgg  versus  v\=R/Zc  is  shown  in  Figure  33  for  the 
same  condition  from  which  Figure  10  is  obtained.  For  the  longer  cables,  energy 
is  deposited  in  the  load  most  rapidly  when  the  cable  is  nearly  matched,  as  has 
been  indicated  by  D.  Tasca,  et  al  (reference  5).  When  a  longer  cable  is  mismatched, 
the  longer  wave  transit  times  induce  lengthly  oscillatory  ringdowns  in  load  currents 
so  that  each  reflected  wave  spends  significant  time  traveling  betwen  loads.  Thus, 
the  time  for  the  load  to  absorb  energy  from  the  cable  increases  rapidly  as 
varies  from  the  matched  condition.  Figure  34  shows  a  plot  of  normalized  t5g  (the 
time  required  for  load  to  absorb  50%  of  the  total  energy)  for  a  family  of  curves 
with  the  dielectric  constant  between  1  and  9  and  for  cable  lengths  between  0.3m 
to  30  meters.  Figure  35  shows  the  normalized  plot  for  a  family  of  values  of  the 
drive  pulse  width. 

For  electrically  short  cables,  the  load  current  tends  to  converge  to  a 
wave  shape  proportional  to  the  driver  pulse 

1L1<t) 

as  indicated  in  Section  2-3.4  (Equation  2.52)  so  that  tgg  convertes  to  T  .  It 
is  interesting  to  note  that  for  the  cable  where  xis  approximately  equal  to  T  , 
tgg  increases  as  R  gets  small.  This  is  exactly  the  point  of  caveat  of  the  deriva¬ 
tion  of  Section  2-3.3,  namely,  that  when  the  cable  delay  is  a  significant  fraction 
of  pulse  duration,  mismatches  produce  significant  ringdowns.  The  plots  of  tgg 
and  t^g  exhibit  rather  broad  based  minimums  in  Figures  33  to  35.  These  minima 
are  a  function  of  driver  pulse  shape,  as  well  as  the  line  parameters  t  andY)  . 
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2-4.2  Square  Pulse  Drivers  on  a  Matched  Line 

Equations  2.35  and  2.36  show  that,  for  a  resistively  loaded  lossless  trans¬ 
mission  line,  the  load  currents  consist  of  sums  of  fundamental  incident  pulses  which 
are  reflected  and  re-reflected  off  the  mismatches.  When  the  transmission  line  is 
matched,  r0  =  0  and  rs  =  0  so  that  the  load  currents  become  simply: 

i^U)  =  ia(t)  =  ~2~  |^(^>  ‘  QU-Tq  cosb-t)|  (2.59) 

1 ( t)  =  ib(t)  =7~"|cl(t-T0  cose)  -  q(t-x ) |  (2.60) 

These  equations  are  useful  to  study  the  effects  of  angle  of  incidence  9,  line 
delay  t,  and  X-ray  propogation  time  t0.  Important  relationships  between  t  and  t0 
will  emerge  that  help  place  the  concept  of  an  electrically  short  line  on  firmer 
footing. 


For  simplicity,  assume  i o ( t )  to  be  square  pulse  driver,  of  strength  Ip 
and  width  Tp.  Then,  for  an  arbitrary  angle  of  incidence  9,  line  delay  r,  the  load 
current  waveforms  are  provided  in  Figures  36  and  37  by  direct  application  of  Equa¬ 
tions  2.59  and  2.60.  The  peak  load  currents,  rise  times,  full  width-half  maximum 
durations,  and  total  energies  of  each  waveform  are  provided  in  Tables  2  and  3 
where  the  load  resistances  are  equal  to  the  line  characteristics  impedance  because 
of  the  assumed  matched  condition.  Except  when  t=t0  and  9=0  the  waveforms  are 
trapezoidal.  In  addition,  the  rise  time  of  the  load  currents  is,  in  general,  a 
function  of  t,  Tp,  and  (t0  cos  9).  The  distinction  between  what  constitutes  an 
electrically  long  line  and  a  short  line  is  very  clearly  defined*.  Technically 
short  and  long  line  waveforms  are  dependent  on  the  load  location,  since  the  load 
current  rise  times  and  pulse  amplitudes  are  a  function  of  where  the  load  is  with 
respect  to  illumination.  Load  1  (z=o)  is  assumed  to  be  illumined  before  Load  2 
(z=s)  as  in  the  derivations  of  Section  2-3.  However,  the  line  can  always  be  con¬ 
sidered  short  (when  the  line  is  matched)  for  a  square  pulse  if 

Tp  >  t  +  tq  cose 

*~fVrThfs~special  case,  see  Section  2-3.3  for  the  general  case 
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Figure  36.  Load  current  waveforms  for  a  square  pulse  driver; 
matched  line  of  delay  i  and  lengths 
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Figure  37.  Load  currents  for  a  square  pulse  driver 
on  a  matched  lossless  air  line  of  delay 
1  length  S  and  ••=() 
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Table  3.  Incident  pulse  parameters  for  a  square  pulse  driver  air  line  with 
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and  long  if 


TP  -  cose 

Intermediate  pulse  widths,  namely 

t-tq  cose  <  Tp  <  t+tq  cose 

result  in  either  "long"  or  "short"  line  responses,  depending  on  load  location.  The 
total  incident  charge  on  one  of  the  two  loads  may  be  noted  as 

% =  s  W2 

and  the  peak  load  current  is  directly  proportional  to  Qq  in  all  cases  when  t<t0. 

The  peak  current  is  also  a  function  of  line  wave  velocity  vp  when  the  line  is  long, 
and  the  faster  the  wave  velocity,  the  higher  the  peak  current,  since  the  accreted 
charge  tends  to  arrive  faster.  When  t  <  tq,  the  ratio  of  short  line  peak  current  to 
long  line  peak  current  in  both  loads  is  somewhat  less  than  the  ratio  of  line  lengths, 
because  the  longer  lines  tend  to  keep  somewhat  more  charge  in  the  lines  as  a  func¬ 
tion  of  time. 

In  the  limit,  when  vQ  =  vp  and  o  =  0,  Table  3  shows  that  the  peak  line 
current  is  independent  of  x  in  load  2  and  is  exactly  Tp  in  duration  regardless  of 
line  length.  In  load  1,  the  peak  of  the  long  line  pulse  is  ID  Tp  vp/4  where  vp  =  vQ. 
Since  Tp  <  2x  =  2S/v0,  the  peak  current  is  less  than  or  equal  to  SIp/2  in  load  1  for 
the  long  line. 

Thus,  it  is  seen  that  for  very  long  lines,  significant  changes  in  pulse 
amplitude  and  duration  can  be  expected  to  occur  when  the  line  is  not  instantaneously 
illuminated,  i.e.,  when  e  f  -n/2.  In  particular  load  1,  first  illuminated,  receives 
more  energy  than  load  2  when  both  loads  are  matched.  However,  this  is  compensated 
somewhat  in  that  the  pulse  duration  in  load  1  is  longer  than  in  load  2. 

Because  of  the  complexity  of  the  line  voltage  equation  (Equation  2.13), 
little  has  been  said  about  voltage  other  than  to  recommend  a  numerical  Fourier 
inverse  for  its  computation.  However,  the  results  of  Figures  36  and  37  provide  some 
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interesting  insights  to  the  terminal  open  circuited  voltage  for  a  square  pulse 
driver.  With  Qq  =  S*D  Tp  =  the  final  accreted  charge,  the  open  circuit  cable  voltage 
is  ultimately 


w  QpZC  SID  TP 

OC  T  T  ZC 


(2.61) 


On  the  first  reflection,  the  open  circuited  line  terminal  voltage  is  simply 

Voc  =  2ia(t)ZC  at  Load  1  (2.62a) 

Voc  =  21b(t)zc  at  Load  2  (2.62b) 


where  the  superscript  denotes  the  voltage  due  to  the  first  incident  pulse.  For  both 
loads , 


V' 


oc 


oc 


except  in  the  case  of  load  2  for  a  long  line  (Figure  36f).  With  x  >  (Tp+-t0  cose) 


peak  V*  =  2 
oc 


iini. 


Sir,  T 


2(t-i0  cose) 


zc>. 


D  'P 


Z  =  V 
c  oc 


(2.63) 


That  is  to  say  that  the  peak  open  circuited  first  reflection  terminal  voltage  at 
load  2  is  greater  than  the  voltage  which  results  from  Equation  2.61. 

Because  of  reflections,  this  voltage  peak  will  occur  at  terminal  1  on  the 
second  pulse.  Since  the  cable  is  assumed  lossless,  the  cable  will  ring  ad  infinitum 
and  never  settle  to  a  quiescent  Voc. 

It  is  stated  without  proof,  that  when  e  =  tt/2,  Voc  =  Qg  Zc/x  is  the  maxi¬ 
mum  peak  terminal  voltage  and  also  the  peak  line  voltage. 
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SECTION  3 

GENERAL  MULTIWIRE  TRANSMISSION  LINE  THEORY 


3-1  BASIC  EQUATIONS 

A  multiwire  transmission  line  is  defined  as  a  collection  of  parallel 
conductors,  all  of  which  are  separated  from  a  reference  conductor  by  dielectric 
material.  For  the  cases  of  interest  in  this  report  the  transmission-line  is  as¬ 
sumed  to  be  perturbed  by  a  current  source  term  driven  by  incident  ionizing  radia¬ 
tion.  In  the  following  sections  we  will  sketch  the  origin  of  the  defining  equations 
for  a  multiwire  transmission  line  and  in  Appendix  I,  II  and  III,  we  will  consider 
detailed  solutions  and  the  meaning  of  the  various  terms  that  appear  in  these 
equations. 

Let  us  consider  a  length  aX  of  the  multiwire  transmission-line  consisting 
of  N  conductors  and  a  reference  conductor.  Each  conductor  will  have  an  inductance 
and  a  capacitance  with  respect  to  the  reference  conductor  and  with  respect  to  every 
other  conductor,  which  we  denote  by  LAX  and  CaX,  where  L  and  C  are  inductance  and 
capacitance  per  unit  lenqth  respectively.  The  voltage  drop  across  AX  is  equal  to 
the  product  of  the  inductance  of  this  length  of  the  transmission  line  and  the  time 
rate  of  change  of  current  summed  over  all  conductors  and  the  reference  conductor: 


A  V .  * 

1 


•(LiJ);X 


(3.1) 


The  decrease  in  current  across  this  length  of  line  is  equal  to  the  current  that  is 
shunted  across  the  distributed  capacitance  plus  the  current  source.  The  change  in 
the  current  is  given  by  the  capacitance  multiplied  by  the  time  rate  of  change  of 
the  voltage  plus  the  current  source: 
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where  ( I q) i  is  the  current  source  per  unit  length  for  wire  i. 
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In  the  limit  aX-<0,  these  equations  become 
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The  starting  point  of  the  previous  derivation  was  the  assumption  that  the 
line  coul d  be  represented  by  a  collection  of  segments  each  with  a  distributed  capa¬ 
citance  and  inductance  along  the  line.  A  more  rigorous  derivation  would  start  with 
Maxwell's  equations  appropriate  tc  guided  waves  and  by  making  the  assumption  that 
the  electric  and  magnetic  fields  are  both  perpendicular  to  the  direction  of  propo- 
gation  and  to  each  other  (TEM  mode)  the  previous  set  of  differential  equations 
would  be  obtained. 

3-2  NATURE  OF  THREAT 

In  this  section  we  derive  the  features  of  multiwire  transmission  line  solu¬ 
tions  which  are  necessary  for  arriving  at  a  practical  engineering  method  of  specify¬ 
ing  current  drive  levels  for  the  multipoint,  multiwire  current  injection  technique 
which  has  been  developed.  As  has  been  noted  in  the  Section  3-1,  the  threat  to  the 
interface  circuits  is  the  energy  arriving  at  the  terminations  due  to  the  distributed 
current  drivers  llong  the  caDle  length.  The  exact  wave  shape  at  the  termination 
depends  not  only  on  the  parameters  of  the  cable  (current  drivers,  characteristic 
impedance  matrix,  propoqation  velocity),  but  also  on  the  termination  impedances 
themselves.  Similarly,  the  actual  energy  which  is  delivered  to  any  particular  load 
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also  depends  both  on  the  cable  parameters  and  the  terminations.  Thus,  to  ascertain 
the  magnitude  of  the  threat  to  any  particular  interface,  it  is  necessary  to  consider 
how  the  signal  arriving  at  that  interface  depends  on  all  the  other  parameters  of  the 
problem.  That  is,  a  worst  case  approach  must  be  developed  for  specifying  this  threat. 
While  it  is  relatively  easy  to  find  upper  bounds  for  this  threat  t put  all  current 
into  a  single  termination) ,  this  naive  approach  can  lead  to  excessive  levels  which 
will  be  above  any  physically  realistic  possibility.  Instead,  it  has  been  decided  to 
look  at  the  multiwire  interface  threat  in  detail  and  extract  a  physically  realizable 
worst  case;  that  is,  not  just  an  upper  bound,  but  a  least  upper  bound. 

It  is  the  purpose  of  this  section  to  explore  the  nature  of  interface 
threats  due  to  multiconductor  transmission  lines  driven  by  X-ray  generated  current 
drivers.  The  physically  realizable  worst  case  will  be  determined.  An  engineering 
prescription  will  be  developed  for  specifying  the  current  drive  levels  required  to 
realize  this  worst  case  when  implementing  multipoint,  multiwire  current  injection 
on  component  boxes. 

The  procedures  for  developing  the  prescription  relies  on  the  method  of 
successive  complication.  An  ideal  baseline  model  is  chosen  for  developing  the  neces¬ 
sary  concepts.  A  definitive  solution  is  presented  for  this  model.  Departures  from 
this  ideal  are  then  discussed  by  successive  relaxation  of  constraints.  The  same 
procedure  is  followed  within  the  baseline  model.  First,  the  nature  of  the  interface 
threat  for  a  single-wire  to  shield  transmission  line  is  reviewed.  Next,  the  simplest 
multiwire  line  which  includes  coupling  is  discussed,  the  two-wire  to  shield  line. 
Finally,  the  general  multiwire  line  is  discussed,  and  the  results  related  to  the  pre¬ 
vious  discussion.  The  engineering  prescription  then  becomes  apparent. 

3-2.1  The  Basic  Model 

The  model  chosen  for  the  baseline  investigation  is  that  of  an  ideal, 
loss-less  m  .1 ticonductor  transmission  line  which  may  be  described  by  the  equa¬ 
tions  of  "quasi-TEM1  transmission  line  theory.  Thus,  the  assumptions  are: 

i)  The  propogation  may  be  described  by  Equations  3.5  and  3.6. 

ii)  The  line  is  combed  and  uni  form, i .e . ,  C  and  C  are  independent 
of  X  and  t. 

iii)  The  current  driver  In  is  uniform  along  the  line.  This  cor¬ 
responds  to  photon  illumination  uniform  and  perpendicular 
to  the  cable  run. 


iv)  All  modes  of  the  line  have  the  same  propagation  velocity  v 
(see  Appendix  I).  This  assumption  simplifies  the  analysis 
and  is  not  far  from  reality.  Mathematically,  this  require¬ 
ment  places  a  constraint  on  the  transmission  line  matrices, 
L  and  C,  i .e. , 


”“  =  C  L  =  \  JL  (3.7) 

V 

where  ||  is  the  unit  matrix. 

v)  For  simplicity,  we  assume  the  line  to  be  symmetrically  ter¬ 
minated.  That  is,  the  reflection  matrices  for  both  ends 
are  identical. 

vi )  Only  passive  resistive  terminations  are  considered.  Re¬ 
active  terminations  dissipate  no  energy  and  would  increase 
the  scope  of  the  investigation  significantly. 

vii)  All  terminations  are  wire-to-shield  only.  No  cross-wire 

terminations  are  considered.  This  corresponds  to  the  situa¬ 
tion  of  primary  interest  in  the  realistic  deployment  of  such 
cables. 

With  the  above  assumptions,  a  closed  form  analytic  solution  of  the  problem 
may  be  obtained.  The  details  of  this  solution  may  be  found  in  Appendix  II,  while 
details  of  the  energy  deposited  in  the  loads  may  be  found  in  Appendix  III.  These 
solutions  will  be  used  extensively  in  the  following  discussion.  Additionally, 
generalizations  of  the  solution  will  be  used  when  discussing  the  relaxation  of  some 
of  the  basic  model  constraints. 

3-2.2  Single  Wire  Line 

The  general  solutions  given  in  Appendix  II  for  a  multiwire  line  of  course 
also  gives  the  single  wire  (to  shield)  solution.  Let  R  be  the  termination  resis¬ 
tance,  fc  the  characteristic  impedance  of  the  line,  and  n  =  R/Zc  the  impedance 
ratio.  The  reflection  coefficient  r  is  given  by 
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(3.8) 


and  varies  from  -1  for  a  shorted  termination,  to  +1  for  an  open  termination.  For 
a  current  driver  Ip  (amps/m)  in  an  electrically  short  cable  the  current  I  in  the 
load  is  given  by 

Kt)  =  -  2S-  O-r) 

£  v  2n  ( aQ( t-2nr )  +  faQ  (t-(2n+l)x) 

n=0 
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where  t  =  s/v  is  the  cable  propagation  time,  v  is  the  propagation  velocity  in  the 
cable,  and  s  is  the  cable  length.  The  basic  pulse  ia  is  given  by 


AQ(t)  =  q(t)  -  q(t-i) 


(3.10) 


where 


ID(t -)  dt  " 


(3.11) 


is  the  basic  charge  transfer  per  unit  length  developed  by  the  current  driver  I[)(t). 
This  is  the  same  relation  developed  in  Section  2-3.3  for  the  special  case  of  an 
electrically  shot  cable  with  ;'o  =  i's- 


For  the  purposes  of  further  discussion,  it  is  assumed  that  the  time  history 
of  the  current  driver  Iq  is  a  square  pulse.  Evaluation  of  the  expressions  for  other 
pulse  shapes  will  be  obvious.  Let  Tp  be  the  width  of  the  square  pulse,  and  Iq  the 
amplitude  of  the  pulse,  as  shown  in  Figure  30a  .  The  charge  transfer  q(t)  ramps 
up  to  a  total  transfer  Q0= I pTp  in  a  time  Tp,  as  shown  in  Figure  3o|.  .  The  basic 
pulse  in  the  termination  aQ  is  somewhat  different  depending  on  whether  the  cable 
is  long  compared  to  the  pulse  (x-Tp).  For  a  long  cable,  aQ  ramps  to  a  peak  value 
Q0  =  I pTp  in  a  time  Tp,  remains  at  this  value  until  r  and  then  ramps  back  to  zero 
at  a  time  Tp  +  x.  For  a  short  cable,  AQ  ramps  to  the  smaller  value  of  IpT  in  a  time 
i,  remains  at  this  level  until  Tp,  and  ramps  back  to  zero  at  a  time  Tp  +  i.  These 
wave  forms  may  be  seen  in  Figure  38c  and  38d  . 

The  important  point  is  that  the  waveform  aQ  (t-m),  which  appears  in 
Equation  3.9,  is  non-zero  only  in  the  interval  nTp<t^(n+l)  Tp+t.  Outside  of  this 
interval  A.Q(t-nr)  is  zero.  Thus,  the  successive  terms  in  the  series  of  Equation 
3.9  appear  at  later  and  later  times.  Thus,  the  terms  are  easily  interpreted  as 
the  successive  reflections  from  the  terminations ,  with  the  primary  reflection  occur¬ 
ring  at  either  end,  followed  by  propagation  to  the  opposite  end  and  reflection  back 
to  the  termination  of  interest.  As  each  pass  through,  the  termination  dissipates 


sone  of  the  energy  of  the  pulse;  the  successive  terms  arrive  with  smaller  and  smaller 
amplitude  (i.e.,  the  pulse  is  multiplied  by  a  power  of  p). 

The  primary  pulse  arriving  at  the  termination  is  most  easily  seen  by  con¬ 
sidering  the  special  case  of  matched  termination  (n=l»  r=0)  and  t  is  small  for  any  value 
of  n,  then  the  infinite  series  of  Equation  3.9  reduces  to 


I(t)  =  --§*  'Q(t)  (3.12) 

which  is  the  same  as  derived  in  Section  2-3.3  with  o-n / 2.  Thus,  the  expression 
reproduces  the  well-known  fact  that  matching  terminations  give  rise  to  no  reflections. 
The  current  which  flows  in  the  termination  is  directly  obtained  from  the  information 
about  AQ  previously  presented.  It  takes  on  its  maximum  value  only  for  i>Tp.  For 
mismatched  terminations  (T70),  the  current  waveform  takes  on  a  different  character 
depending  on  whether  there  is  a  high  mismatch  (n>l,  r<0)  or  a  low  mismatch  (n<l,  r>0) . 
The  first  case  has  the  character  of  an  over-damped  ring  down  of  a  charged  capacitor, 
while  the  second  case  waveform  has  the  character  of  an  under-damped  R,L,C  circuit. 
These  features  of  the  termination  waveform  may  be  seen  in  Figures  39  and  40  which 
give  the  current  waveforms  for  a  moderately  long  cable  (i=Tp)  for  significant  high 
mismatch  (R=10Zo)  and  low  mismatch  (R=0.1Zo). 

Several  other  features  of  the  solution  should  be  noted.  The  load  current 
scales  directly  with  the  source  current  10.  It  also  scales  with  the  phase  velocity 
v.  For  long  cables,  the  current  scales  as  vQ0=vI0Tp,  while  for  short  cables,  it 
scales  as  vIqi-Iqs.  The  load  current  depends  on  the  load  resistance  only  through 
the  dimensionless  ratio  n=R/Z0.  Finally,  the  time  dependence  of  the  solution  in 
terms  of  a  dimensionless  time  variable  t  =t/i  depends  only  on  the  pulse  length 
through  the  dimensionless  ratio  c=Tp/i.  Taken  together,  this  means  that  the  solu¬ 
tion  space  depends  only  on  the  parameters  n  (equivalently  r)  and  r.,  and  that  other 
solutions  may  be  obtained  by  scaling.  Of  course,  the  load  voltage  may  be  obtained 
from  the  current  by  multiplication  by  R. 

While  load  current  and  voltage  are  of  some  concern  for  semi-conductor 
terminations,  it  is  recognized  that  many  major  failure  mechanisms  are  associated 
with  the  total  energy  deposited  in  the  load,  and  the  rate  at  which  the  energy  is 
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Figure  40.  Terminal  response  -  low  mismatch 
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deposited.  Thus,  while  the  current  in  the  load  reaches  a  maximum  as  the  length 
of  the  cable  is  increased  IDTp  for  a  matched  load),  the  length  of  time  the  cur¬ 
rent  flows  continues  to  increase  as  the  cable  length  increases.  This  implies  that 
longer  cables  will  deliver  more  energy  to  the  load  even  though  the  peak  current 
and  voltage  may  be  no  worse  than  for  a  shorter  cable.  The  development  of  worst 
case  considerations  for  semi-conductor  loads  demands  that  these  energy  concerns  be 
addressed.  A  central  feature  of  the  approach  to  .  •*  discussed  below  is  a  careful 
accounting  for  worst  case  energy  capability.  This  criterion  is  taken  to  be  of  prime 
importance.  The  features  of  this  energy  capability  are  now  discussed  for  the  single¬ 
wire  transmission  line.  An  understanding  of  this  capability  will  be  central  to  the 
approach  for  multi -wire  lines. 


In  Appendix  III,  a  derivation  has  been  given  for  the  energy  which  is 
deposited  in  the  loads  on  a  multiwire  transmission  line.  For  a  single  wire  long 
line  ( x>Tp )  this  expression  reduces  to 


where 


(3.13) 


(3.14) 


and  C  =  capacitance  per  unit  length.  In  these  equations,  the  symbols  have  the 
same  meaning  as  above,  i.e.,  c  =  Tp/t  is  the  characteristic  time  ratio,  Q0  >  I[)Tp 
is  the  charge  transfer  per  unit  length,  C  is  the  capacitance  per  unit  length  of 
the  line,  s  is  the  line  length,  and  r  the  reflection  coefficient.  Equation  3.14 
may  also  be  written  in  terms  of  the  characteristic  impedance  of  the  line  Zc  =  v^C"^, 
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For  shorter  lines,  (x  <  Tp)  the  expressions  for  the  energy  deposited 
become  more  complicated,  as  is  noted  in  Appendix  III.  However,  the  resulting  expres¬ 
sion  always  gives  rise  to  an  energy  less  than  that  implied  by  Equation  3.13  ,  i.e., 
more  and  more  positive  terms  are  subtracted  from  this  basic  quantity.  Furthermore, 
as  the  termination  impedence  becomes  very  large  (R-*»,r-»-  + 1 ) ,  all  of  the  expressions 
tend  to  Em  .  no  matter  what  the  ratio  of  T  to  t.  For  further  discussion,  it  is 

ITlaX  p 

therefore  sufficient  to  consider  the  case  Tp  <  x,  as  given  be  Equation  3.13  . 

The  maximum  energy  available  Emax  is  simply  the  amount  of  energy  developed  by  charging 
up  the  capacitance  of  the  cable  with  the  total  transfer  Q  .  Naturally,  this  quantity 
scales  linearly  with  the  length  of  the  cable.  It  is  thus  possible  to  discuss  the 
per  unit  length  energy  capability  of  the  cable.  Note  that  this  quantity  depends  not 
only  on  the  magnitude  of  the  current  driver  IQ  (through  Qq  =  I^Tp)  but  also  on  the 
cable  capacitance  (or  equivalently  the  characteristic  impedence  Z  ).  Note  that  for 
very  long  cables  U<<1),  the  energy  delivered  to  the  load  is  nearly  independent  of 
the  load  resistance.  Shown  in  Figure  41  is  the  normalized  energy  E/Emax  as  a  function 
of  n  =  R/Zq  for  various  values  of  c.  Note  that  this  quantity  is  a  monotically  increas¬ 
ing  function  of  the  load  resistance  R.  Shown  in  Figure  4?  is  the  same  information 
displayed  in  absolute  units  of  E  which  shows  the  increase  in  energy  with  cable  length. 

Not  only  is  the  total  energy  of  importance  for  semi-conductor  considera¬ 
tions,  but  also  the  energy-time  history.  Shown  in  Figure  43  are  the  energy  time 
histories  for  three  identical  cables  with  low,  matched,  and  high  impedence  termina¬ 
tions  (n=0.1,  1,  10).  The  solution  is  for  a  long  cable  (x  =  Tp),  so  the  final  energies 
are  very  nearly  the  same.  As  is  evident  from  the  figure,  the  energy  is  deposited  in 
the  load  much  more  rapidly  for  the  matched  load  as  compared  to  the  mismatched  loads. 
While  no  convenient  analytic  expression  has  been  developed  for  the  length  of  time 
required  to  deposit  a  fixed  fraction  of  the  energy,  it  is  evident  from  the  Equation 
3.9  that  the  matched  condition  clearly  gives  rise  to  peak  power  in  the  load. 

Indeed,  for  very  long  lines,  since  the  final  deposited  energy  is  nearly  independent 
of  the  load,  it  is  clear  that  the  matched  load  represents  the  worst  case.  Indeed, 
the  time  to  deposit  the  energy  in  the  load  is  minimized  by  matching  the  line.  Shown 
in  Figure  44  is  a  summary  of  the  time  to  deposit  one  half  the  final  energy  as  a 
function  of  termination  impedence.  The  deep  minima  for  most  line  lengths  is  apparent. 
This  plot  was  developed  from  an  extensive  numerical  evaluation  of  the  solution  for 
the  indicated  parameter  range. 
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Figure  42.  Energy  vs.  R/Z 


Sufficient  information  has  now  been  presented  about  single  wire  solutions 
to  proceed  to  the  simplest  multi-wire  model,  that  of  the  two-wire-to-shield  trans¬ 
mission  line. 

3-3  TWO  WIRE  ANALYSIS 

It  is  of  interest  to  consider,  in  a  multi-wire  bundle,  how  one  can  couple 
current  from  one  wire  to  another.  The  obvious  place  to  study  such  effects  is  the 
two-wire  bundle.  In  the  Appendix  we  have  considered  the  analytic  analysis  of  the 
two-wire  bundle  and  have  shown  how  the  two-wire  bundle  reduces  to  single  wire  solu¬ 
tions.  Such  a  reduction  is  extremely  useful  in  terms  of  conceptualizing  what  is 
happening  in  the  coupling,  but  the  obvious  way  to  demonstrate  the  coupling  is  to 
present  the  results  of  the  current  and  energy,  for  example,  that  appear  in  the  two 
wires  when  only  one  is  driven.  The  current  that  flows  in  the  un-driven  wire  is  all 
coupled  current.  We  have  analyzed  various  two-wire  bundles  in  order  to  consider 
various  coupling  effects. 

The  analytic  solution  that  is  presented  in  Appendix  II  assumes  that  there 
is  only  one  propagation  velocity.  In  general,  in  a  multiwire  bundle,  there  can  be 
many  propagation  velocities,  which  are  a  source  of  coupling  from  one  wire  to  another. 
While  it  is  beyond  the  scope  of  this  report  to  detail  the  consequences  of  this  coupl¬ 
ing,  some  calculation  have  been  done  on  a  two-wire  bundle  with  two  different  phase 
velocities.  The  analysis  seems  to  indicate  that,  particularly  at  low  impedance  ter¬ 
mination,  there  can  be  a  significant  coupling  effect  due  to  the  different  phase 
velocities.  The  remainder  of  the  calculations  have  been  performed  on  two-wire  bundles 
with  one  phase  velocity.  One  two-wire  bundle  was  constructed  using  the  same  capaci¬ 
tance  matrix  as  the  multi -phase  velocity  bundle  and  by  choosing  the  inductance  matrix 
in  such  a  way  as  to  give  one  phase  velocity.  The  other  two-wire  bundle  was  constructed 
in  such  a  way  as  to  represent  a  seven-wire  bundle  where  six  of  the  seven  wires  are 
tied  together.  Let  i  represent  the  wire  of  interest  and 


Vt  -  V  if  i 

Vi  ®  V  i  =  l  (3.16) 
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where  Vi  is  the  voltage  of  the  itn  wire  with  respect  to  the  shield.  The  relationship 
between  charge  and  voltage  is  again 


C.  -V  . 
ij  J 


If  we  sum  over  all  i,  i  t  fc  we  have 


(3.17) 


=  Yc.  .V  . 

^  1 3  3 

i  ,  j  At  * 
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and 


N 


■  Zct/j 
j=l  J  J 


if  we  think  in  terms  of  two  wires,  Q-j  and  Q2  and  let 


(3.18) 


(3.19) 
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(3.20) 
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we  have 


V  C11V1  +  C12V2 


Q2*  C12  V1  +  C22  V2 


(3.21) 


where 


CH=  Cll 


C12  SCilj  ~  C21  "2Cji, 

vfcc,  (3-22i 

22  lT^t 1J 

We  have  considered  a  seven  wire  bundle  and  picked  out  the  wire  with  the  largest 
driver  so  that  i=3  in  the  seven  wire  bundle. 


In  Figure  45  and  46  we  plot  the  ratio  of  E/E for  the  bundles  TW2  and 
TW3  respectively.  In  this  ratio,  E  is  the  energy  dissipated  across  each  termina¬ 
tion  for  the  un-driven  wire  and  EMAX  is  the  maximum  energy  capability  across  each 
termination  for  the  driven  wire.  The  quantity  E/EMAX  is  plotted  as  a  function  of 
R2  for  various  values  of  ,  where  R2  is  the  ratio  of  the  termination  impedance 
to  the  diagonal  characteristic  impedance  of  the  un-driven  wire  and  R1  is  the  same 

ratio  for  the  driven  wire.  These  plots  indicate  that  the  coupled  energy  from  the 

2 

driven  to  the  un-driven  wire  approach  a  maximum  value  of  C  times  the  maximum 
value  of  the  energy  of  the  driven  wire,  C  is  the  coupling  factor  defined  in  Appendix 
V.  It  is  important  to  realize  that  the  quantity  C  is  always  much  less  than  one 
for  any  realizable  bundle. 


3-4  SEVEN  WIRE  ANALYSIS 

The  seven  wire  line  analysis  has  been  performed  in  order  to  demonstrate 
the  basic  tenets  of  the  previous  sections  when  applied  to  a  larger  bundle  of  wires. 

We  have  considered  various  distributions  of  seven  wire  drivers  which  are  typical  of 
groups  of  seven  wires  in  a  larger  bundle  (37  wire  bundle).  In  these  calculations 
we  have  singled  out  the  wire  in  that  distribution  with  the  largest  driver  and  varied 
the  terminating  resistance  of  that  wire  as  a  function  of  the  terminating  resistances 
of  the  other  six  wires.  In  Appendix  VII  we  present  the  results  of  these  calculations 
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One  basic  tenet  of  our  analysis  is  that  one  can  think  of  each  wire  in  a 
multiwire  bundle  as  a  single-wire  transmission-line  driven  by  an  effective  drive 
Ip  given  by 

N 

iD  ■  It  +  EciiIii 

(3.23) 

Where  I  j  is  the  driver  for  the  wire  of  interest,  wire  *  ,  I,.  are  the  drivers  for 
all  the  other  wires,  and  are  the  coupling  terms. 


C. 

1 


=zi  */<ZiiZu>5 


(3.24) 


The  charcteristic  impedance  matrix  is  given  by  Z  and  N  is  the  number  of  wires  in 
the  bundle.  We  can  compute  the  total  energy  dissipated  for  a  single-wire  trans¬ 
mission  line,  Emax,  driven  by  this  effective  driver  using  the  results  of  Appendix 
III  and  compare  with  the  actual  energy  dissipated  from  the  numerical  calculations. 
Shown  in  Figure  47-49  are  the  results  of  these  calculations.  The  labelling  conven¬ 
tion  is  detailed  in  Appendix  III.  On  these  plots,  is  the  maximum  energy  capa¬ 

bility  per  termination  for  the  wire  with  the  largest  driver  (including  the  coupling), 
R1  is  the  ratio  of  the  termination  impedance  to  the  diagonal  element  of  the  character¬ 
istic  impedance  for  the  wire  with  the  largest  driver,  and  is  the  same  ratio  for 
all  the  other  wires.  It  is  clear  from  these  plots  that  EMAX  represents  the  upper 
limit  for  the  energy. 


3-5  LOSSY  CABLE  TRANSMISSION  LINE 

It  is  of  interest  to  consider  dissipative  effects  in  evaluating  the 
current  that  flows  through  the  terminations  for  a  cable  driven  by  X-rays.  In  general, 
for  any  real  cable  there  will  be  a  resistance  per  unit  length  along  the  length  of 
the  cable  due  to  skin  depth  losses  and  a  conductance  per  unit  length  across  the  di¬ 
electric  due  to  the  dark  conductivity  of  the  dielectric.  In  the  following  para¬ 
graphs,  we  will  discuss  how  to  incorporate  frequency  dependent  conductance  and 
resistance  into  the  transmission-line  equations  and  present  our  solutions  of  these 
equations. 
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Figure  47.  E(JOUEES)  vs. 


MAX 


Figure  49.  E(JOULES)  vs. 


If  we  allow  our  caDle  to  be  lossy,  the  following  transmissionline  equa¬ 
tions  are  appropriate: 


.-CSV 

at 


£JL  =_,-£4  -  GV  +  I, 


(3.25) 


For  convenience  we  will  limit  ourselves  to  single-wire  cables  and  so  drop  all 
vector  notations.  In  these  equations  R  is  the  resistance  per  unit  length  and  G 
is  the  conductance  per  unit  length.  The  method  of  solution  is  the  same  as  indi¬ 
cated  in  Appendix  II,  and  so  we  arrive  the  following  fourier-transformed  equations. 


-  iwC?  - 
||  =  iwcfy  +  1 


(3.26) 


At  this  point,  we  will  allow  R  and  G  to  be  functions  of  frequency,  i.e. 


R  =  R(W) 
G  -  G(W) 


(3.27) 


Using  the  same  assumptions  as  made  in  Appendix  II  we  arrive  at  the  following 
equations: 


d2fr  =  -k2^  +  I iWL-R  L 

J  o  l  D 


(3.28) 


where 


2  w 

k  =  — 2  -  RG  +  iw(RC  +  GL). 
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(3.29) 
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These  equations  can  be  solved,  as  in  Appendix  II,  for  I(w,X)  and  V(w,X).  In  par- 

% 

ticular,  the  solution  for  1(0, w)  is  as  follows 


a,  w) 

f(o,w)  =  ■  D 
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(2.30) 


The  notation  in  this  solution  is  the  same  as  has  been  used  previously  with  the 
following  expression  for  ZQ,  the  characteristic  impedance 

Zq  =  k/(wc  +  iG)  (3.31) 

Since  k  and  ZQ  are  functions  of  w,  it  is  not  possible  to  obtain  directly  I(o,t), 
given  by  the  following  equation 

-foo 

I(0,t)  "  2¥  jdw  l(o,w)e"lwt  (3.33) 

—  oo 

We  have  performed  the  required  integrations  numerically  using  a  fourier-transform 
code  and  will  present  the  results  of  this  analysis  in  Appendix  VIII. 
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SECTION  4 
MODEL  VALIDATION 


4-1  INTRODUCTION 

The  key  to  successfully  developing  the  desired  threat  definition  is  an 
accurate  analytical  tool  which  can  predict  not  only  photon  responses  (i.e.,  the  de¬ 
sired  invariant  current  source  terms)  but  also  the  cesired  electrical  parameters 
required  for  the  maximum  energy  analysis.  Existing  tools,  such  as  the  CHIC  code, 
have  been  developed  to  provide  such  information.  An  interesting  characteristic  of 
these  programs,  which  is  associated  with  the  physics  of  the  direct  drive  phenomenon, 
is  that  the  calculated  electrical  parameter  and  current  driver  definitions  are  derived 
from  the  same  set  of  physical  and  material  parameters  contained  in  the  subject  cable. 
The  fact  that  these  entities  are  not  mutually  independent  of  one  another  provides  a 
fortuitous  situation  which  can  be  exploited  to  advantage. 

Since  the  successful  development  of  the  maximum  energy  concept  rests  with 
the  accuracy  of  such  tools  as  the  CHIC  code,  it  is  highly  desirable  to  obtain  an 
independent  validation  of  the  accuracy  of  such  codes  before  committing  extensive 
effort  based  on  their  utilization.  Such  a  validation  is  described  in  this  section. 

The  interrelationship  of  the  current  driver  and  electrical  parameter  definitions 
provides  an  excellent  double  check  on  the  code  validation  conclusions.  The  cable 
type  which  will  be  selected  is  also  the  baseline  multiwire  cable  type  which  has  been 
used  as  the  model  for  the  development  of  the  multiwire  analysis  technique.  An  excel¬ 
lent  choise  for  a  candiate  cable  type  is  the  AWG  24  (Raychem  44/0411  space  wire),  37 
wire  foil/film  (Copper/Mylar)  shielded  cable  identical  to  that  used  in  the  DSCS  III 
Spacecraft  and  in  the  DSCS  III  Harness  Photon  Test. 

4-2  CABLE  PARAMETER  SENSITIVITY-DRIVERS 

The  very  large  parameter  space  associated  with  any  possible  multiwire  cable 
configuration  requires  a  considerable  reduction  if  any  meaningful  predictive  capa¬ 
bility  is  to  be  realized.  That  is  to  say,  we  would  like  to  isolate  those  parameters 
which  predominantly  affect  the  response  of  a  multiwire  cable  configuration  to  inci- 
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dent  ionizing  radiation.  In  order  to  effect  this  reduction,  we  have  employed  the 
CHIC  computer  code  to  perform  a  large  number  of  calculations  in  various  multiwire 
cable  configurations.  The  CHIC  code,  details  of  which  have  been  presented  elsewhere, 
(8-10) ^ s  uniqueiy  suited  to  perform  this  analysis  since  it  is  capable  of  handling  an 
arbitrary  cross  sectional  geometry  and  material  configuration. 

For  this  particular  cable  configuration,  the  parameter  space  that  is  per¬ 
tinent  to  an  SGEMP  analysis  would  include  the  following  variables: 


Individual  Wires 

- 

Composition,  Plating,  Size 

Wire  Insulation 

- 

Composition,  Thickness,  Dielectric 
Constants,  Density 

Liner 

- 

Composition,  Thickness,  Density 

Overshiel d 

- 

Composition,  Plating,  Thickness 

Geometric 

Shield  configuration  relative  to 
wire  bundle,  relative  positions  of 
wires,  total  number  of  wires  in 
bundle. 

Our  approach  in  this  sensitivity  study  is  to  consider  only  those  variables 
which  we  feel  might  realistically  be  considered  in  any  mul ti -conductor  bundle.  In 
addition,  we  have  considered  how  one  might  reduce  these  variables  even  further  for 
the  purpose  of  a  sensitivity  study.  One  example  of  this  is  the  assessment  problem 
associated  with  multi-layered  dielectrics  surrounding  the  individual  wires  of  the 
bundle,  with  frequency  dependent  dielectric  constants.  We  have  investigated  the 
possibility  of  representing  such  a  configuration  with  one  dielectric  with  an  effec¬ 
tive  dielectric  constant.  The  outcome  of  this  investigation^^  has  indicated  that 
it  is  perfectly  adequate,  for  the  purposes  of  a  sensitivity  study,  to  model  a  multi¬ 
layered  dielectric  by  a  single  dielectric  with  an  effective  dielectric  constant. 

Such  a  modeling  will  certainly  indicate  the  trends  in  the  response  as  a  function 
of  the  various  parameters.  On  the  other  hand,  if  one  wants  to  pinpoint  the  response 
of  a  given  multiwire  bundle,  one  has  to  take  into  account  the  details  of  the  multi¬ 
layered  dielectric  configuration.  Such  an  analysis  has  been  performed  and  it  has 

been  found  that  a  modeling  of  the  double  layered  dielectric  is  critical  to  accurately 

(121 

calculating  the  response  of  the  multiwire  bundle. 
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4-2.1  Sensitivity  Study 


The  following  parameters  have  been  considered  in  our  sensitivity  study: 

•  Wire  Size  -  16-26  AWG 

•  Wire  Plating  -  Sn,  Ni 

•  Shield  Composition  -  Cu,  A1 

•  Liner  Thickness  -  1,2  mil 

•  Angle  of  Incidence  of  Radiation 

•  Geometric  Configuration 

•  Shield  Configuration 

4-2. 1.1  Geometric  and  Shield  Configuration 

One  question  of  interest  is  the  dependence  of  the  maximum  response  in  a 
bundle  to  the  number  of  wires  in  the  bundle.  Shown  in  Figure  50  is  a  comparison  of 
a  61-wire  bundle  and  a  7-wire  bundle.  The  two  bundles  are  identical,  except  for  the 
number  of  wires.  Since  the  largest  response  is  always  on  the  wire  nearest  to  the 
shield,  we  do  not  expect  there  to  be  a  significant  difference  in  the  maximum  response 
from  a  61-wire  bundle  to  a  7-wire  bundle.  In  fact,  our  calculations  indicate  that 
the  maximum  response  for  the  61-wire  bundle  is  slightly  smaller  than  the  7-wire 
bundle.  For  this  reason,  we  have  considered  the  7-wire  bundle  to  be  a  canonical 
configuration  representing  any  close-packed  mul ti -conductor  bundle  (see  Figure  51 
and  have  performed  the  majority  of  our  computations  on  it.  We  have  also  considered 
various  packing  arrangements  of  the  wires  in  the  bundle;  the  7-wire  bundle  representing 
hexagonal  close-pack,  and  have  found  very  little  difference.  An  example  of  a  hexa¬ 
gonal  close-pack  and  circular  pack  bundle  of  wires  is  shown  in  Figure  52.  In  addition, 
we  have  found  very  little  difference  in  the  response  between  a  circular  and  hexagonal 
shield,  and  so  have  performed  most  of  our  calculations  on  a  circular  shield  configura¬ 
tion.  An  example  of  a  hexagonal  shield  is  shown  in  Figure  53. 

4-2. 1.2  Wire  Size 

We  have  considered  four  different  wire  radii,  ranging  in  size  from  8.0 
to  25.4  mils.  Shown  in  Figure  54  are  worst  case  response  versus  wire  gauge.  Since 
a  positive  response  indicates  that  the  dominant  effect  is  wire  emission,  we  expect 
that  the  bundle  with  larger  radii  wires  will  have  the  larger  responses,  and  this 
is  verified  by  the  calculations. 
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Figure  50.  Comparison  of  61-wire  and  7-wire  bundle 
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Figure  51.  Geometric  model  used  for  parametric  study 
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Figure  52.  Comparison  of  hexagonal  and  circular  pack 
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Figure  54.  Worst  case  response  versus  wire  gauge  -  2  mil  liner 


We  have  considered  Sn  and  Ni  wire  plating  and  Cu  and  A1  shield  composi¬ 
tion.  Since,  again,  the  dominant  effect  is  wire  emission,  we  expect  the  largest 
response  to  occur  for  that  configuration  which  has  the  largest  Z  wire  plating  and 
the  smallest  Z  shield  composition  and  this  is  verified  by  the  calculations,  as 
indicated  in  Figure  54. 

4-2 .1.4  Liner  Thickness 

The  effect  of  the  thickness  of  the  dielectric  liner,  which  is  bonded  to 
the  overshield,  on  the  response  of  the  wires  in  the  bundle  has  been  considered. 

A  thinner  liner  would  result  in  the  shield  emission  being  more  effective,  thereby 
reducing  the  worst  case  response  for  a  given  bundle.  We  have  considered  one  and 
two  mil  liners.  Shown  in  Figure  55  is  a  typical  distribution  of  drivers  for  these 
two  liners.  It  is  clear  that  the  response  is  indeed  smaller  when  the  liner  is  half 
as  thick,  but  not  significantly  smaller. 

4-2.1  .5  Angle  of  .Incidence  of  Radiation 

Shown  in  Figure  56  are  the  distribution  of  drivers  for  various  angles  of 
incidence.  An  angle  of  incidence  of  0°  gives  rise  to  the  worst  case  driver. 

4-2.2  Conclusions 

For  the  purposes  of  comparing  one  bundle  configuration  to  another,  we 
will  use  as  the  criterion  the  maximum  response  in  a  given  bundle.  For  the  range 
of  parameters  considered  in  this  study,  the  bundle  configuration  that  resulted  in 
the  largest  maximum  response  had  the  following  parameters: 

(a)  16  AWG  wire  size, 

(b)  Sn-plated  wire, 

(c)  Al  shield, 

(d)  2  mil  liner  thickness,  and 

(e)  0°  angle  of  incidence. 

The  maximum  response  for  this  bundle  was  computed  to  be  6.12  x  lO"1^ 

2 

(coul/cm)/(cal/cm  ).  If  we  take  this  configuration  to  be  the  baseline  configura¬ 
tion,  the  following  specific  conclusions  can  be  drawn  from  the  sensitivity  study: 
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Figure  55.  7-wire  response 


Figure  56.  7-wire  response 
HO’10  coul/cm/cal/cm2] 
effect  of  incidence  angle 


(1)  The  maximum  response  is  a  monotonically  increasing  function 
of  wire  radius  for  all  material  configurations. 

(2)  Replacing  an  aluminum  shield  with  a  copper  shield  gives  at 
least  a  factor  of  0.64  reduction  in  driver  from  the  baseline 
response . 

(3)  Replacing  the  tinned  wires  by  Ni-plated  wires,  gives  a  factor 
of  0.58  reduction  in  driver  from  the  baseline  response. 

(4)  Replacing  the  tinned  wires  by Ni -plated  wires,  and  replacing 
an  aluminum  shield  by  a  Cu  shield  gives  a  factor  of  0.26 
reduction  in  driver  from  the  baseline  response. 

(5)  For  all  material  configurations,  a  small  liner  gives  a 
slightly  decreased  response  for  the  worst  case  driver. 

(6)  The  angle  of  incidence  of  the  radiation  has  only  a  small  effect 
on  the  worst  case  driver,  and  0°  incident  is  the  worst  case. 

These  results  indicate  the  relative  insensitivity  of  a  mul ti -conductor 
bundle  to  realistic  variations  in  the  parameter  space.  In  general,  if  we  take  the 
maximum  response  as  an  indicator,  the  maximum  variation  is  less  than  a  factor  of 
ten.  Although  this  analysis  has  mainly  been  carried  out  for  one  typical  blackbody 
spectrum,  our  calculations  have  also  indicated  that  the  response  increases  with 
blackbody  temperature,  a  result  consistent  with  the  absence  of  gaps  between  all 
conductor-di electric  interfaces . 

4-3  MODEL  DRIVERS 

This  section  will  present  the  results  of  a  complete  analysis  of  the  elec¬ 
trical  parameters  and  current  drivers  associated  with  7  and  37  foil/film  multiwire 
cables.  This  analysis  has  been  performed  using  the  CHIC  code,  det.ils  of  which 

(o\ 

have  been  considered  elsewhere/  1  The  CHIC  code  has  the  capability,  not  only 
of  calculating  the  current  sources  associated  with  Photon  Direct  Drive  SGEMP,  but 
also  all  the  parameters  which  electrically  identify  a  given  multiwire  cable.  The 
electrical  parameters  can  be  conveniently  represented  in  matrix  form  and  are  completely 
analogous  to  the  capacitance  and  inductance  of  a  coaxial  cable. 

4-3.1  Foil/Film  Multiwire  Cable  -  Specification 

We  will  be  considering  a  7  and  37  wire  bundle  in  this  section.  These  two 
Lundies  are  identical,  except  for  the  difference  in  the  number  of  wires.  Each  multi¬ 
wire  cable  consists  of  N  individually  unshielded  wires  within  a  single  foil/film 
overshield.  In  this  section,  we  will  specify  completely  the  geometrical  and  material 
configuration  which  has  been  used  in  these  calculations. 
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4-3. 1.1  Individual  Wires 

The  individual  wires  are  Raychem  AWG  #24  space  wire  which  consists  of 
pure  tin  coated,  copper  wire  covered  by  a  combination  cross-linked,  extruded  poly- 
alkene  and  cross-linked,  extruded  polyvinyl  dene  fluoride  dielectric  (Kynar).  The 
extrusion  process  used  for  the  Raychem  space  wire  is  such  that  there  are  no  gaps 
between  the  dielectric  and  the  conductor.  Shown  in  Figure  57  is  a  model  of  the 
Raychem  space  wire  with  appropriate  labels.  Table  4  is  a  table  of  the  correspond¬ 
ing  parameters  associated  with  the  Raychem  space  wire  which  serves  as  inputs  to 
the  CHIC  code. 

4-3. 1.2  Foil/Film  Overshield 

The  foil/film  overshield  consists  of  a  Sun  Chemical  #1132  film-foil  lami¬ 
nate  made  of  an  electro  deposited  Copper  foil  on  a  polyester  (Mylar)  film.  Figure 
58  is  a  model  of  the  overshield  with  appropriate  labels  and  Table  5  is  a  table 
of  the  associated  parameters  which  also  serve  as  inputs  to  the  CHIC  code. 

4-3 .1.3  Geometric  Model 

We  have  assumed  in  these  calculations  that  the  individual  wires  in  the 
multiwire  bundle  are  hexagonally  close-packed  and  that  the  overshield  is  circular. 
Shown  in  Figures  59  and  60  are  the  geometric  models  assumed  in  the  calculations 
along  with  the  labeling  of  the  individual  wires  in  the  bundle.  This  labeling  con¬ 
vention  will  be  used  throughout  this  section,  both  in  the  definition  of  the  capaci¬ 
tance  and  inductance  matrices  and  in  the  current  drivers. 

4-3.2  Incident  Spectrum 

The  current  drivers  have  been  computed  using  the  P I  MBS - 1 A  spectrum  shown 

in  Figure  61.  We  have  assumed  a  20  mil  Mylar  and  1  mil  A1  window  between  the 

photon  source  and  the  multi -wire  cable,  and  will  present  all  current  drivers  nor- 
2 

malized  to  1  cal/cm  incident  on  the  cable  overshield. 

4-3.3  Electrical  Parameter  Definition  of  a  Multiwire  Bundle 

There  are  two  capacitive  matrices  which  are  relevant  to  multiwire  bundles 
of  N  wires.  One,  which  we  will  call  Maxwell's  capacitance  coefficients,  is  defined 
by  the  following  matrix  equation: ^ 

N 

‘  2>ijVj  .  i-1 . N 
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Figure  57.  RAYCHEM  space  wire  *24  AWG 


Table  4.  Parameters  associated  with  RAYCHEM  space  wire  ~  24  AWG 


MATERIAL 

RADIUS 
(mils ) 

DENSITY 

3 

(gm/cm  ) 

RELATIVE 

DIELECTRIC 

CONSTANT 

A 

(TIN  COATED  COPPER) 

r1  =  13.0 

8.96 

B 

(POLYALKENE) 

r2  =  16.5 

1  .08 

2.25 

C 

(KYNAR) 

r3  =  19.5 

1.73 

7.0 
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8  Figure  58.  Foil /film  overshield 


Table  5.  Parameters  associated  foil/film  overshield 


THICKNESS 

DENSITY 

DIELECTRIC 

MATERIAL 

(MILS) 

(GM/CM3 ) 

CONSTANT 

A  (COPPER  FOIL) 

tx  =  0.7 

8.96 

— 

B  (MYLAR) 

t2  =  °-92 

1  .35 

2.9 

I 

i 
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Figure  59.  7-wire  multiwire  bundle  with  labelling  convention 


Figure  60.  37-wire  multiwire  bundle  with  labelling  convention 
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Intensity  Normal  to  Anode  (per  electron,  per  2*  steradians) 


Figure  61.  Sandyl  calculation  of  PIMBF  spectrum  incident  on  cables 
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where  Qi  is  the  charge  on  the  ith  conductor,  Vj  is  the  voltage  of  the  jth  conductor 

with  respect  to  the  shield,  and  C-.  are  Maxwell’s  capacitance  coefficients.  C.  .  is 

J  th  ^  th 

the  ratio  of  the  charge  induced  on  the  i  conductor  to  the  potential  on  the  j 

conductor  when  all  other  conductors  are  at  zero  potential.  An  alternative  capacitance 

matrix  which  can  be  determined  from  C. .  will  be  labeled  (CM)...  In  this  definition, 

(C^)^.  is  the  ratio  of  the  charge  on  the  i  conductor  to  the  potential  on  the  i 

conductor  when  all  other  conductors  have  no  charge.  (CM) ■  •  (i^j)  is  the  ratio  of 
t  h  ^ 

the  charge  on  the  i  conductor  to  the  potential  difference  ( V ^ - V ^ )  when  all  other 
conductors  have  no  charge. 

The  inductance  matrix  is  determined  by  computing  the  capacitance  matrix 
of  the  identical  multiwire  bundle  without  dielectric.  The  inductance  matrix  L^. 
is  determined  by  the  following  equation: 


where  C  is  the  speed  of  light  in  vacuum. 

2 

The  propagation  matrix  (y  is  determined  once  and  are  known 
and  is  defined  by  the  following  equation: 

N 

(y  ^ij  ”  LikCkj 

k=l 

9  Q 

This  matrix  is,  in  general,  not  diagonal;  but  when  brought  to  diagonal  form  (y  )^; 
the  diagonal  elements  are  related  to  the  multiwire  phase  velocities  by  the 
following  equation: 


A  knowledge  of  L..  and  C-.  can  also  be  used  to  determine  the  characteristic  impe- 
i  J  i  J 

dance  matrix  Z.j  of  a  multiwire  bundle  which  can  best  be  represented  symbolically 

as:  _ _ 

Z  =  V  L/C 

where  2,  L,  and  l  are  the  impedance,  inductance,  and  capacitance  matrices,  respec¬ 
tively. 
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4-3. 3.1  Electrical  Parameter  Results 


In  this  section,  we  present  the  complete  7-wire  results  and  a  limited 
amount  of  37-wire  results.  In  Tables  6-10  are  given  the  complete  set  of  para¬ 
meters  for  the  seven  wire  bundle.  In  Tables  11-15  are  given  the  diagonal  elements 
of  all  the  matrices  for  the  37-wire  bundle. 

4-3.4  Fall-Off  of  Capacitance  Matrix 

For  a  given  wire,  i ,  in  a  multiwire  bundle,  we  do  not  calculate  all  of 
the  elements  of  the  capacitance  matrix  C^,  j=l,  2,  ...,  N,  but  only  those  elements 
which  correspond  to  wires  that  are  nearest  to  wire  i.  This  is  justified  in  light 
of  the  fact  that  C7- .  is  very  small  for  wire  j  far  away  from  wire  i,  and  is,  there¬ 
fore,  for  the  purposes  of  any  calculation  zero.  To  justify  this  assumption,  we 
have  computed  C^j,  j =1 ,  ...,  N,  for  the  37-wire  bundle.  Shown  in  Table  16  are 
the  results  of  this  calculation. 

4-3.5  Current  Driver  Results 

Using  the  geometrical  and  material  configuration  described  in  Section 
4-3.1  and  the  x-ray  spectrum  described  in  Section  4-3.2,  we  have  computed  the  cur¬ 
rent  drivers  for  the  7-  and  37-multiwire  bundles.  The  results  are  indicated  in  Figure 

62  and  Tables  17  and  18.  The  largest  driver  for  the  7-wire  bundle  is  =3  with  a  driver 
-10  2 

of  +4.14  x  10  coul/cm/cal/cm  .  The  largest  negative  driver  for  the  37-wire 

-10  2 

bundle  is  -14  with  a  driver  of  -1.27  x  10  coul/cm/cal/cm  . 

4-3.6  Comparison  to  Experiment 

It  is  of  interest,  once  the  model  drivers  and  associated  electrical 
parameters  have  been  computed,  to  compare  results  with  actual  multiwire  bundles. 

The  particular  quantities  that  have  been  measured  are  the  Norton  equivalent  current 
drivers  for  a  37-wire  bundle  and  various  capacitance,  wire  to  shield  and  wire  to 
wire,  for  that  bundle.  When  comparing  theory  and  experiment,  however,  it  is  impor¬ 
tant  to  keep  in  mind  the  fact  that  the  geometrical  model  assumes  a  close-packed 
configuration  whereas  there  is  a  certain  tendency  for  the  wires  to  move  around  in 
the  actual  bundle.  It  is  difficult  to  make  a  detailed  comparison  wire  by  wire  for 

this  reason.  Our  model  predictions  always  have  the  wire  with  the  largest  driver 

-10  2 

at  the  outer  edge  of  bundle  with  the  peak  \4  X  10  coul/cm/cal/cm  ,  and  the 

-10  2 

experimental  results  give  the  largest  driver  ^6  X  10  coul/cm/cal/cm  .  The  wire 
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Table  6.  Capacitance  matrix  C  (pf/m)  for  7-wire  bundle 


248.08 

-31.86 

-31.86 

-31.86 

-31.86 

-31.86 

-31 .86 

-31.86 

225.76 

-37.83 

-0.55 

-0.08 

-0.55 

-37.83 

-31  .86 

-37.83 

225.76 

-37.83 

-0.55 

-0.08 

-0.55 

-31.86 

-0.55 

-37.83 

225.76 

-37.83 

-0.55 

-0.08 

-31  .86 

-0.08 

-0.55 

-37.83 

225.76 

-37.83 

-0.55 

-31 .86 

-0.55 

-0.08 

-0.55 

-37.83 

225.76 

-37.83 

-31 .86 

-37.83 

-0.55 

-0.08 

1 

-0.55 

-37.83 

225.76 

Table  7.  Capacitance  matrix  CM  (pf/m)  for  7-wire  bundle 


207.19 

130.10 

130.10 

130.10 

130.10 

130.10 

130.10 

130.10 

203.05 

128.63 

110.01 

107.64 

110.01 

128.63 

130.10 

128.63 

203.05 

128.63 

110.01 

107.64 

110.01 

130.10 

110.01 

128.63 

203.05 

128.63 

110.01 

107.64 

130.10 

107.64 

110.01 

128.63 

203.05 

128.63 

110.01 

130.10 

110.01 

107.64 

110.01 

128.63 

203.05 

128.63 

130.10 

128.63 

110.01 

107.64 

110.01 

128.63 

203.05 

Table  8.  Inductance  matrix  L  (nH/m)  for  7-wire  bundle 


187.72 

40.58 

40 . 58 

40.58 

40.53 

40.58 

40.58 

40.5  8 

152.11 

31  .79 

12.71 

10.05 

12.71 

31.79 

40.58 

31  .  79 

152.11 

31  .  79 

12.71 

10.05 

12.71 

40.58 

12.71 

31.79 

152.11 

31.79 

12.71 

10.05 

40.58 

10.05 

12.71 

31.79 

152.11 

31.79 

12.71 

40.58 

12.71 

10.05 

12.71 

31  .79 

152.11 

31 .79 

40.58 

31.79 

1  2  .  71 

10.05 

12.71 

31  .79 

152.11 

23 


,=  2  20  -2 

Table  9.  Propagation  matrix  y  (10  (m/s)  )  for  7-wire  bundle 


3881 .29 

6.25 

6.25 

6.25 

6.25 

6.25 

6.25 

206.46 

3063.77 

-37.46 

-9.88 

-3.27 

-9.88 

-37.46 

206.46 

-37.46 

3063.77 

-37.46 

-9.88 

-3.27 

-9.88 

206.46 

-9.88 

-37.46 

3063.77 

-37.46 

-9.88 

-3.27 

206.46 

-3.27 

-9.88 

-37.46 

3063.77 

-37.46 

-9.88 

206.46 

-9.88 

-3.27 

-9.88 

-37.46 

3063.77 

-37.46 

206.46 

-37.46 

-9.88 

-3.27 

_ 

-9.88 

-37.46 

3063.77 

Table  10. 


Phase  velocities  (e-values  V 


(108  m/s) 


for  7-wire  bundle 


1.61 
1.84 

1 .80 

1.80 
1.79 

1 .81 


24 


C<1 ) 


247.29  pf/m 


C<1,32> 


Table  17.  7-wire  current  drivers 


with  the  largest  driver  tends  to  be  on  the  outside  of  the  bundle  in  such  a  way  that 

there  is  a  separation  between  it  and  the  nearest  neighbors.  We  believe,  therefore, 

that  the  difference  between  the  model  predictions  and  the  experimental  data  are  due 

to  the  fact  that  the  wires  are,  in  fact,  not  totally  close  packed.  The  model  drivers 

range  from  -1.27  X  10”^  to  +4.2  X  10 coul/cm/cal /cm^ .  Another  way  to  compare 

theory  and  experiment  is  to  compare  the  bulk  current  drivers  defined  as  the  algebraic 

.0 

sum  for  the  bundle.  The  model  predicts  a  bulk  current  ID^4.2  X  10  .  The  difference 

p 

between  the  model  predictions  and  the  experimental  data  can  more  vividly  be  seen  by 
the  two  histograms  plotted  in  Figure  63  and  64.  In  the  model,  there  are  more  negative 
drivers  in  the  interior  of  the  bundle  due  to  the  close  packing  of  the  wires.  This  is 
what  gives  rise  to  a  smaller  bulk  current. 

Measurements  have  been  made  of  the  capacitance  matrix  of  a  37-wire  bundle, 
and  it  is  of  interest  to  compare  these  results  with  the  model  predictions.  It  is 
important  to  note  that  whatever  is  the  reason  for  differences  between  the  model  and 
the  data  as  far  as  the  current  drivers  are  concerned  can  also  be  responsible  for 
difference  in  the  capacitance  matrix.  We  can  think  of  the  capacitance  matrix  as  de¬ 
fining  the  geometry  and  dielectric  properties  of  the  bundle.  Table  19  is  a  comparison 
of  the  experimental  data  of  various  elements  of  the  capacitance  matrix  and  the 
corresponding  model  prediction.  In  all  cases  the  data  are  somewhat  smaller  than  the 
model.  The  difference  is  consistent  with  the  fact  that  the  wires  are  not  as  close 
packed  as  the  model  assumes. 
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Figure  63.  37-wire  model  -  histogram  of  normalized  distribution 


Table  19.  Comparison  of  capacitance  matrix  elements  -  model  and  experiment 


WIRE  TO  SHIELD  | 

WIRE 

TO  WIRE 

CAPACITANCE 

CAPACITANCE 

(PF/M) 

MODEL 

EXP.  DATA 

(PF/M) 

MODEL 

EXP.  DATA 

C(1,1) 

168 

90 

C  ( 1 , 2 ) 

112 

86 

C(2,2) 

124 

90 

C ( 2, 3) 

106 

85 

C (3, 3) 

124 

89 

C (3,4) 

106 

85 

C ( 4/ 4) 

124 

90 

C  ( 4 , 5 ) 

106 

85 

C ( 5, 5) 

124 

89 

C  (  5 , 6  ) 

106 

85 

C  ( 6 , 6 ) 

124 

90 

0(6,7) 

106 

86 

C ( 7, 7) 

124 

90 

C  ( 7, 8) 

78 

67 

C ( 8/ 8 ) 

120 

95 

C  ( 8, 9) 

81 

65 

C (9,9) 

119 

95 

C  ( 9 , 1 0 ) 

81 

65 

C(  10, 10) 

119 

96 

C(10,11) 

81 

65 

C(11,11) 

119 

95 

C ( 1 1  ,  1  2) 

81 

65 

C(12,12) 

119 

96 

C  ( 1  2 , 1  3) 

81 

65 

C  ( 1 3,  1  3) 

119 

96 

C( 13,14) 

106 

83 

C  ( 1 4, 1  4) 

125 

91 

C(14,15) 

80 

64 

C(  15,1  5) 

124 

90 

C(15,16) 

80 

63 

C(16,16) 

124 

90 

C( 16,17) 

80 

64 

C  ( 1  7, 1 7) 

124 

91 

C (17,18) 

80 

64 

C  ( 1  8 , 1  8) 

124 

91 

C(18,19) 

80 

64 

CC 19,19) 

124 

96 

C ( 1 9, 20) 

70 

60 

C ( 20 , 20 ) 

124 

108 

C ( 20 , 2 1 ) 

103 

84 

C  (  2 1 , 2  1 ) 

124 

110 

C ( 21 , 22) 

79 

68 

C(22,22) 

125 

107 

C (22,23) 

103 

86 

C  C  23, 23) 

124 

112 

C ( 23, 24) 

79 

68 

C(24,24) 

125 

110 

C ( 24 , 2  5 ) 

103 

86 

C(25,25) 

124 

112 

C ( 2  5 , 26) 

79 

69 

C(26,26) 

125 

108 

C (26,27) 

103 

86 

C(27,27) 

124 

109 

C ( 2  7, 28 ) 

79 

68 

C (28,28) 

125 

108 

C ( 28, 29) 

103 

86 

C(29,29) 

124 

113 

C ( 29, 30) 

79 

69 

C ( 30 , 30 ) 

125 

1 1 1 

C  (  30, 31 ) 

103 

86 

C(31,31) 

125 

1 1  5 

C(  31 , 32) 

1 1  1 

87 

C ( 32 , 32) 

196 

1 1  3 

C (32, 33) 

105 

65 

C ( 33, 33) 

196 

109 

C ( 33, 34) 

105 

63 

C (34, 34) 

196 

113 

C ( 34, 35) 

105 

65 

C (35,35) 

196 

112 

C  (  35 , 36) 

105 

65 

C ( 36, 36) 

196 

1 1  1 

C ( 36 , 3  7) 

105 

65 

C(37,37) 

196 

115 

C  (  34 , 37) 

100 

63 

C (33,36) 

100 

61 

C ( 32, 35) 

100 

62 

C ( 1 , 32 ) 

99 

60 

C( 1,33) 

99 

60 

C  ( 1 , 34  ) 

99 

60 

C ( 1 , 35 ) 

99 

60 

C ( 1 , 36) 

99 

60 

C( 1, 37) 

112 

60 

33 


SECTION  5 

CIT  DESIGN  AND  DEVELOPMENT 


5-1  INTRODUCTION 

The  ground  characteristic  of  the  photon  direct  drive  threat  has  been  shown 
in  Sections  2  and  3  to  be  one  of  a  large  differential  mode,  wi re-to-shield  transmission 
line  type  electrical  surge.  This  type  of  threat  produces  a  unique  energy-time,  trans¬ 
mission  line  controlled  waveshape  at  electronic  interfaces  and,  it  is  to  this  threat 
profile  that  satellite  electronics  are  hardened.  Hence,  the  primary  requirement 
imposed  on  any  candidate  C.I.T.  technique  is  that  it  provides  a  meaningful  SGEMP 
hardness  verification  test  with  respect  to  correctly  simulating  the  electrical  energy¬ 
time  nature  of  the  SGEMP  threat.  There  are,  however,  a  number  of  additional  design 
constraints  which  are  imposed  at  the  satellite  system  program  level  which  impact  the 
specific  method  and  configuration  that  can  be  used  to  meet  the  C.I.T.  specification. 
These  constraints  are  both  of  a  practical  nature  as  well  as  those  associated  with 
actual  spacecraft  hardware.  Specifically,  the  C.I.T.  must  be  compatible  with  the 
SGEMP  hardening  designs  employed  in  actual  satellite  systems,  it  must  be  of  a  low  risk 
approach  which  provides  individual,  positive  control  of  currents  delivered  to  each  pin, 
and  finally,  be  practical  for  user  test  evaluation. 

An  additional  requirement  is  associated  with  the  correct  electrical  parameter 
definition  of  the  maximum  SGEMP  energy  threat  which  exists  in  a  cable  bundle,  since 
it  is  these  parameters  which  must  be  used  in  a  meaningful  hardness  verification  test. 
Furthermore,  the  maximum  energy  capability  of  a  particular  wire  is  not  only  controlled 
by  the  magnitude  of  the  wire's  current  source  term  but  is  also  controlled  by  the 
wire-to-shield  characteristic  impedance. 

The  situation  is  further  complicated  by  the  fact  that  since  satellite 
cables  are  comprised  of  many  unshielded  wires  enclosed  by  a  common  bundle  shield, 
significant  interwire  coupling  occurs  which  alters  the  equivalent  electrical  para¬ 
meters  of  the  transmission  line  type  threat. 

From  a  threat  specification  standpoint,  one  could  adopt  an  approach  based 
on  defininq  the  maximum  threat  posed  to  the  component  box  pins  in  terms  of  the 
total  cable  core  current  and  some  nominal  or  worst  case  characteri stic  impedance. 
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This  is  an  extemely  conservative  approach  which  is  unrealistic  but  mathematically 
representative  of  an  upper  limit.  This  approach,  however,  is  completely  intolerable 
from  a  satellite  system  standpoint  since  is  represents  a  significant  hardware  impact 
design  solution  to  providing  adequate  hardening  against  a  realistic  threat  level. 

5-2  DESIGN  CONSIDERATIONS 

Since  the  SGEMP  threat  in  present  day,  design  hardened,  satellite  system  is 
associated  with  photon  direct  drive  effects  in  electrical  cables,  the  foundation  for 
the  overall  hardening  designs  for  component  boxes  is  based  upon  the  energy-time, 
transmission  line  nature  of  the  direct  drive  threat.  However,  since  system  weight, 
power,  cost  and  reliability  are  of  utmost  concern,  these  hardening  designs,  by  neces¬ 
sity,  are  of  a  hardware  efficient  nature  for  both  new  and  existing  design  component 
boxes.  Thus,  SGEMP  hardening  is  implemented  on  a  "specifically  where  required"  basis 
rather  than  legislated  on  an  "across  the  board"  basis. 

Since  the  severity  of  the  SGEMP  threat  is  controlled  by  energy-time  trans¬ 
mission  line  characteristics  and  specific  wire  length,  hardening  requirements  are 
imposed  based  on  actual  cable  lengths  in  the  satellite  system  for  each  interface  cir¬ 
cuit  as  well  as  the  quantitatively  defineable  electrical  loads,  not  only  at  the  circuit 
of  concern,  but  at  the  interfacing  circuit  at  the  opposite  end  of  the  wire.  In  many 
cases,  for  example,  the  presence  of  a  particular  low  impedance  electrical  load  or  the 
addition  of  a  single  zener  diode  at  one  end  of  a  short-to-medium  length  wire  is  suf¬ 
ficient  to  provide  protection  to  the  interfacing  circuits  at  both  ends  of  the  line. 

For  longer  length  lines,  SGEMP  hardening  is  implemented  at  both  ends  of  the  line 
consistent  with  the  line  length  of  the  transmission  line  type  threat.  Hence,  many 
interface  circuits  in  a  component  box  are  only  SGEMP  hardened  when  all  their  respec¬ 
tive  interface  component  boxes  are  electrically  connected. 

An  additional  design  consideration  is  associated  with  the  physical  nature 
of  the  electrical  systems  (cables  and  connectors)  in  actual  satellites.  Satellite 
cables  are  unlike  re-entry  vehicle  cables  in  that  they  are  not  generally  straight, 
two  connector  cables.  Instead,  satellite  cables  are  typically  complex  configurations 
comprised  of  large  bundles  of  upwards  of  hundreds  of  unshielded  wires  contained  within 
a  common  bundle  shield  with  numerous  branches  and  connector  breakouts,  with  high 
density  (large  number  of  pins)  connectors  used  for  component  box  interfaces.  Further¬ 
more,  the  wiring  in  these  cables  consists  of  wires  with  short  total  lengths  (low 


threat  level  wires  which  require  minimal  SGEMP  hardening  of  their  respective  inter¬ 
face  circuits)  in  a  common  bundle  with  wires  of  long  total  lengths  (high  threat 
level  wires  which  require  extensive  interface  circuit  hardening). 

The  consequence  of  these  design  aspects  is  that  an  SGEMP  simulation  tech¬ 
nique  must  provide  high  simulation  fidelity  of  the  energy-time  transmission  line 
nature  of  the  SGEMP  threat  for  multipin,  simultaneous  drive  of  component  boxes.  It 
also  must  be  capable  of  correctly  driving  a  component  box  whose  SGEMP  protection  is 
contained  in  numerous  other  component  boxes  which  are  not  in  the  test  setup  during 
individual  component  box  testing.  Furthermore,  the  technique  must  be  capable  of 
simul taneously  driving  some  pins  at  high  levels  and  others  at  low  levels  to  accom¬ 
modate  various  distributions  of  long  and  short  wire  threats  at  the  same  connector 
interface.  This  latter  requirement  also  has  a  significant  impact  on  the  risk  potential 
associated  with  any  potential  simulation  technique. 

Since  the  SGEMP  hardness  verification  test  is  typically  performed  in  line 
at  the  latter  portion  of  the  total  qualification  cycle  (generally  before  thermal- 
vacuum  qualification),  a  low  risk  simulation  approach  is  required  to  avoid  inadver¬ 
tent  overstress  levels  which  could  cause  spacecraft  equipment  failures  and  impact 
system  integration  and,  possibly,  flight  schedules.  In  view  of  the  typical  SGEMP 
hardening  designs  employed  in  satellite  systems,  one  cannot  generally  afford  the 
luxury  of  overdriving  some  pins  at  higher  levels  ’n  order  to  assure  that  the  minimum 
required  test  levels  are  properly  established  in  all  pins.  As  such,  the  simulation 
method  must  also  be  capable  of  providing  individual  positive  control  of  each  of  the 
currents  delivered  to  each  of  the  interface  pins  to  eliminate  high  current  (or  voltage) 
hazards  to  both  the  equipment  as  well  as  the  operating  personnel.  An  obvious  addition¬ 
al  requirement,  which  is  more  practically  oriented  for  user  test  evaluation,  is  that 
the  simulation  methods  be  capable  of  being  configured  into  portable,  in-house  hard¬ 
ware  readily  adaptable  to  correctly  interfacing  (both  electrically  and  mechanically) 
with  actual  satellite  component  boxes  and  their  associated  ground  test  equipments. 

5-3  CIT  DESIGN 

The  detailed  design  requirements  were  based  on  the  design  constraints 
described  in  the  previous  section.  In  order  to  achieve  the  desired  electrical  simula¬ 
tion  the  design  and  development  of  the  test  techniques  were  accomplished  in  the 
following  five  subsystems  elements. 
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1)  A  central  source  term  generator  to  reproduce  the  desired 
invariant  current  driver  waveshape  to  drive  a  practical 
multiport  current  source. 

2)  A  set  of  cable  delays  from  the  central  generator  to  the 
multiport  current  source  to  provide  time-phased  photon 
source  term  simulation. 

3)  Discrete  current  source  injection  into  a  controlled  charac¬ 
teristic  impedance  (Z0)  test  cable  for  each  pin  under  test 
to  provide  simulation  of  the  distributed  source  term  nature 
of  the  threat. 

4)  Multiple  cable  line  drive  capability  in  order  to  simulta¬ 
neously  drive  the  number  of  pins  which  corresponds  to  the 
maximum  number  in  the  cable  bundle  which  could  realistically 
be  at  maximum  threat  conditions. 

5)  Provide  a  range  of  termination  impedances  and  drive  levels 
for  the  test  injection  cable  to  simulate  those  corresponding 
to  the  realistic  maximum  threat  when  considering  realistic 
maximum  multiwire  coupling  effects  in  the  actual  multiwire 
cable  which  is  connected  to  the  component  box  being  tested. 

The  central  direct  drive  source  term  generator  provides  a  large  amplitude 
current  pulse  which  has  the  same  time  waveshape  as  the  actual  photon  pulse.  This 
is  the  unit  which  reproduces  the  basic  invariant  current  driver  waveshape  and  is  the 
heart  of  the  C.I.T.  technique.  The  multiport  currents  source  and  time  phasing  cables 
provide  the  correct  distribution  of  the  invariant  current  driver  waveshape  to  the 
pins  of  the  component  box  to  complete  the  simulation. 

A  significant  element  in  the  development  of  the  C.I.T.  approach  was  the 
proper  technique  to  distribute  the  source  term  generator  to  a  multiwire,  controlled 
impedance  network  for  correct  simulation.  This  distribution  depends  on  the  specific 
wire  placement  in  each  cable  manufactured.  Hence,  each  pin  on  the  component  box  must 
be  evaluated  for  maximum  threat  in  order  to  ensure  that  the  complete  interface  is 
survivable  no  matter  what  manufacturing  variations  occur  in  the  production  satellites. 
This  evaluation  must  be  performed  keeping  in  mind  that  the  hardened  level  of  each  pin 
corresponds  to  the  actual  wire  length  in  the  system  and  that  the  indiscriminate  in¬ 
jection  of  current  into  a  high  characteristic  impedance  wire  could  produce  an  over¬ 
stress  condition. 

The  approach  for  distributing  the  source  terms  whi ch  was  more  preferred  from 
a  practical  and  simulation  aspect  was  to  simultaneously  test  the  desired  number  of 
signal  lines  at  maximum  threat  levels  using  controlled  characteristic  impedance 


values  and  current  amplitude  values  corresponding  to  maximum  threat  conditions.  Hence, 
each  pin  should  be  driven  by  its  own  controlled  impedance  pulse  generator  interface 
at  maximum  wire-to-shield  threat  levels  when  considering  realistic  worst  case  multi¬ 
wire  coupling  effects.  Multiwire  transmission  line  analyses  is  required  in  this 
approach  for  the  proper  definition  of  such  effects  and  for  the  development  of  user 
oriented  simplified  techniques  to  readily  define  realistic  test  parameters  to  use 
for  maximum  threat  testing  of  any  specific  cable  design  and  configuration. 

A  simplified  diagram  of  the  multipin  C.I.T.  method  and  equipment  which  was 
developed  is  shown  in  Figure  65.  This  method  is  essentially  a  group  of  "N"  transmis¬ 
sion  lines  driven  by  "S"  current  sources.  Line-to-line  isolation  diodes  are  provided 
as  shown  to  prevent  interwire  bias  loading  and  to  prevent  loading  of  the  injection 
test  cable  by  the  multiport  current  source.  A  bipolar,  central  source  term  generator 
is  formed  by  the  voltage  pulse  across  the  two  output  resistors.  The  multiport  current 
source  is  formed  by  the  drive  resistors  in  series  with  the  isolation  diodes  which  are 
driven  by  the  voltage  pulse  across  the  output  resistors  of  the  central  source  term 
generator.  Photon  time  phasing  simulation  is  accomplished,  quite  simply,  by  driving 
each  section  of  the  line  through  different  lengths  of  coaxial  cable. 

The  pulser  consists  of  the  current  pulse  generator,  various  load  impedances 
and  test  cables.  The  basic  schematic  for  the  C.I.T.  Pulser  is  shown  in  Figure  66. 

A  test  cable  impedance  of  36  ohms  was  used  as  a  result  of  an  analysis  which  indicated 
that  the  average  wire-to-shield  characteristic  impedance  on  typical  DSCS  III  cables 
was  36  ohms.  The  remaining  variables,  such  as,  cable  length  and  termination  impe¬ 
dance  was  specified  by  the  particular  component,  or  connector  under  test.  Simulation 
of  the  distributed  current  source  was  achieved  by  using  the  pulser  to  inject  the  36 
ohm  test  cable  once  per  meter.  A  typical  test  setup  shown  in  Figure  67  used  test 
cable  lengths  of  1/2  meter  between  the  I/F  adaptor  of  the  box  under  test  and  the 
pulser,  1  meter  between  current  injection  points  and  1/2  meter  between  the  last  cur¬ 
rent  injection  point  and  the  termination  load.  A  Tektronix  CT-2  current  probe  was 
used  to  measure  the  current  delivered  into  the  box  under  test.  When  opening  all 
outputs  of  the  pulse  generator  are  either  connected  to  a  test  cable  or  connected  to 
a  short. 

Typical  design  capability  of  this  C.I.T.  equipment  is  shown  below  when 
using  a  10  amp  per  meter  drive  level. 


Bl  POLAR 
CENTRAL 
SOURCE  TERM 
GENERATOR 


Pulse  Generator  Charge: 
Pulse  Shaping  Resistor: 
1  Driver  Max: 


20  kV 
10  ft 

10  Amp/Meter 


Cable  Impedance: 

50  ft 

Cable  Length: 
(T/tp) 

Absolute  Maximum 
Wires  Which  can 

Number  of 
be  Driven 

0.17 

120 

0.5 

60 

1.7 

40 

Figure  68  shows  the  overall  layout  of  the  C.I.T.  equipment  which  was 
designed.  Figure  69  shows  an  overall  view  of  the  attached  pulser  hardware  that  was 
subsequently  developed  based  upon  this  design.  This  hardware  was  then  used  to  test 
actual  spacecraft  boxes  to  demonstrate  the  application  and  utility  of  the  test  tech¬ 
niques  and  pulser  hardware.  The  results  of  these  tests  are  given  in  the  Final 
Technical  Report,  Task  2  of  this  Program  (Reference  13). 

Typical  pulse  current  waveshapes  resulting  from  the  simulated  SGEMP  re¬ 
sponse  of  numerous  cable  lengths  and  which  are  delivered  to  a  large  range  of  compo¬ 
nent  box  pin  impedances  by  the  C.I.T.  method  are  described  here  are  shown  in  Figure  70. 
The  component  box  pin  impedances  "R"  are  expressed  in  terms  of  the  impedance  relative 
to  the  wire-to-shield  characteristic  impedance  (Zc)  of  the  cable  response  being  simu¬ 
lated.  The  cable  electrical  length  threat  t  being  simulated  is  expressed  in  terms 
of  the  cable  propagation  time  compared  to  the  full  width-half  maximum  (FWHM)  of  the 
photon  pulse  being  simulated.  The  time-per-division  "t[)"  and  the  amperes-per-di vision 
"Id"  of  the  oscilloscope  traces  are  expressed  in  time  relative  to  one  another.  Each 
of  the  amplitudes  and  waveshapes  shown  are  in  excellent  agreement  with  those  which 
are  predicted  by  analytical  codes. 
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DESCRIPTION 

HVDC  POWER  SUPPLY 
LVDC  POWER  SUPPLY 
CONTROL  PANEL 
FRONT  PANEL 
LOAD  BOX 
DRIVE  BOX 


Figure  68.  Pulser  layout 
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Figure  69.  Overall  front  view  of  the  C.I.T.  equipment 


SECTION  6 
CONCLUSION 


The  need  for  this  analysis,  design  and  development  effort  is  rooted  in  the 
requirement  for  a  safe  and  reliable  simulation  test  method  to  verify  satellite  elec¬ 
trical  system  operation  in  a  nuclear  X-ray  environment.  A  number  of  fundamental 
transmission  line  response  equations  have  been  developed,  progressing  from  the  simple 
two-wire  case  to  the  multiwire  cable,  which  provide  useful  tools  for  analyzing  cable 
response  in  radiation,  including  the  case  whe^c  the  radiation  propagates  along  the 
length  of  the  cable.  Numerical  and  graphic  results  are  provided. 

For  the.  mul tiwi re  bundle  the  current  driver  description  has  been  reduced 
to  an  effective  single  wire  transmission  line  with  an  effective  current  driver. 

Using  this  effective  current  driver  the  worst  case  current  and  energy  into  an  inter¬ 
face  is  developed  as  a  function  of  line  length  and  termination  impedance. 

The  current  injection  test  method  which  has  been  designed  and  described 
here  represents  what  is  felt  to  be  an  optimum  approach  for  performing  a  meaningful 
multipin  C.I.T.  on  electronic  component  boxes.  The  photon  direct  drive  simulation 
has  been  accomplished  within  the  overall  satellite  program  level  constraints  of 
actual  electronic  equipment. 

The  major  objectives  of  developing  simplified  analytical  and  experimental 
techniques  for  simulation  of  the  photon  direct  drive  threat  to  satellite  electronics 
has  been  met  in  three  areas. 

1.  The  photon/electrical  response  of  multiwire  cables  has  been  determined 
with  respect  to  their  SGEMP  direct  drive  sources,  transmission  line  parameters  and 
realistic  multiwire  coupling. 

2.  The  electrical  test  requirements  for  component  box  current  injection 
t . s  ting  have  been  defined  with  respect  to  the  SGEMP  electrical  energy  distribution 
in  multiwire  cables  and  injection  current  level  and  characteristic  impedance  of 
this  threat. 
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3.  Direct  drive  current  injection  test  simulation  has  been  performed 
using  the  designed  C.I.T.  pulser  on  electrically  active  component  boxes  to  achieve 
a  proper  simulation  of  threat,  at  low  risk,  while  being  compatible  with  actual  hard¬ 
ened  satellite  electronics  and  practical  for  user  test  performances. 
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Appendix  I 


PROPERTIES  OF  THE  TRANSMISSION-LINE  MATRICES 


There  are  two  matrices  which  uniquely  determine 
the  electrical  properties  of  a  multiwire  bundle.  These 
matrices  are  the  capacitance  matrix  C  and  the  inductance 
matrix  L.  The  elements  of  the  capacitance  matrix  C, 
also  called  Maxwell's  capacitance  coefficients,  are 
defined  by  the  following  matrix  equation 

«i  -LCifj  t1-1’ 

j 

where  Qi  is  the  charge  on  the  i'th  conductor,  Vj  is  the 
voltage  of  the  j'th  conductor  with  respect  to  the  shield 
and  Cjj  are  Maxwell's  capacitance  coefficients.  Cjj  is 
the  ratio  of  the  charge  induced  on  the  i'th  conductor  to 
the  potential  on  the  j'th  conductor  when  all  other  con¬ 
ductors  are  at  zero  potential.  The  capacitance  matrix 
has  the  following  properties: 

a)  C  is  symmetric 

b)  The  diagonal  elements  are  all  greater  than  or 
equal  to  zero. 

c)  The  off-diagonal  elements  are  all  less  than  or 
equal  to  zero. 

N 

d)  C.  .  >  0  for  every  i. 

j  =  l  J 


These  properties  can  all  be  proven  from  the  defining 
relationship,  Equation  (1-1),  and  from  the  basic  prin¬ 
ciples  of  electrostatics. 

The  inductance  matrix  L  is  defined  in  a  way 
analogous  to  the  capacitance  matrix 


I,  AX 

J 


(1-2) 


In  this  equation  Ljj  is  the  inductance  coefficient,  i^j 
is  the  magnetic  flux  linking  the  i'th  conductor  and  the 
shield  when  there  is  a  current  Ij  flowing  through  the 
j 1 th  conductor  and  no  current  flowing  through  any  other 
conductor,  and  AX  is  a  short  section  of  length.  This 
equation  merely  states  that  the  inductance  coefficient 
Ljj  is  equal  to  the  flux  linking  the  i’th  conductor 
per  unit  length  of  cable  per  unit  current  flowing  in 
the  j'th  conductor.  The  inductance  matrix  has  the 
following  properties: 


a)  L  is  symmetric. 

b)  All  elements  are  greater  than  or  equal  to 
zero . 

The  propagation  matrix  y2  is  defined  by  the 
following  equation 


~2 

Y 


C  •  L 


(1-3) 


This  matrix  is,  in  general,  not  diagonal;  but  when 
brought  to  the  diagonal  form  by  the  appropriate 
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transformation  matrix,  the  diagonal  elements  are  related 
to  the  multiwire  phase  velocities.  If  we  denote  the 
diagonal  form  of  the  propagation  matrix  as  y^  and  the 
multiwire  phase  velocities  as  V,  then  we  have 


Although,  in  general,  the  multiwire  phase  velocities  can 
all  be  different,  there  are  important  circumstances  in 
which  they  are  all  the  same.  In  particular,  if  there  is 
a  uniform  dielectric  surrounding  all  the  wires  in  the 
multiwire  bundle,  there  is  only  one  phase  velocity  asso¬ 
ciated  with  this  bundle  and  it  will  be  found  that  the 
propagation  matrix  y^  is  diagonal  with  only  one  diagonal 
element  y^.  The  phase  velocity  V  for  this  bundle  is, 
therefore , 


and  we  have 

LC  =  £  L  =  -4  I  (1-6) 

where  I  is  the  identity  matrix. 

Since  the  inductance  matrix  does  not  depend  on 
the  dielectric  properties  of  the  multiwire  cable,  but 
only  on  the  relative  positions  of  conductors,  L  is  the 
same  for  a  given  multiwire  cable  with  or  without  any 
dielectric.  If  we  know  the  capacitance  matrix  C°  of 
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our  multiwire  cable  without  dielectric,  we  can  use  the 
fact  that  there  is  only  one  phase  velocity  for  this 
cable  V=C,  the  speed  of  light,  and  the  relationship 

L  C  *  =  C"5"  L  =  —  T 
C2 


L  = 


(1-7) 


to  determine  the  inductance  matrix  of  our  multiwire 
cable . 


The  characteristic  impedance  matrix  ZQ  of  a 
multiwire  cable  is  determined  once  L  and  C  are  known 
by  the  loxlowing  relation 


ZQ  =  (C'1  L)1/2 


=  C  *(C  L)1/2  =  C  1(y2)1/2 


■«  1  /  r* 


(I»8) 


2 

In  general,  if  y  is  not  a  diagonal  matrix,  but  can  be 
brought  to  the  diagonal  form  by  the  transformation 
matrix  T,  we  have 


XI- 9) 


and,  therefore, 


(1-10) 
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For  the  special  case  when  there  is  only  one  phase 
velocity  V,  we  have 


(1-11) 


RELATION  TO  MEASURABLE  QUANTITIES 

From  the  results  of  the  preceding  section,  it 
is  clear  that  a  complete  description  of  the  electrical 
properties  of  a  multiwire  bundle  is  obtained  once  the 
capacitance  matrix  C  and  is  computed.  As  part  of 
the  computer  code  CHIC,  these  matrices  are  calculated 
and  all  other  relevant  parameters  are  computed.  It  is 
important  to  note  that  a  knowledge  of  C  and  C^  entails 
complete  information  about  the  electrical  properties 
of  a  multiwire  bundle,  and,  therefore,  any  particular 
measurement  made  of  capacitances  and  inductances  can 
be  compared  with  the  theoretical  results  by  appropriate 
manipulation . 

It  is  usually  more  convenient,  when  making 
a  capacitance  measurement,  to  relate  charge  to  the 
potential  difference  between  two  conductors.  If  one 
performs  such  an  exercise  on  a  multiwire  bundle  one 
obtains  a  capacitance  matrix  Cy  which  is  different 
from  the  one  considered  in  the  previous  sections.  In 
the  following  paragraph  we  will  indicate  how  one  obtains 
CM  by  appropriate  manipulation  of  C. 
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T 


The  defining  relation  for  C  is  again 


3  =  C  V 


(1-12) 


We  can  invert  this  relation  and  solve  for  V 


V 


(1-13) 


If  we  specify  the  potential  on  wire  i  as  V ,  let  the 
charge  of  wire  i  be  Q  and  at  the  same  time  let  all  of 
the  other  wires  float  so  that  Qj=0  for  i^j ,  we  have 


V  =  (c~1)iiQ 


(1-14) 


We  can  write  this  relation  in  the  following  way 


Q  = 


-  «*>!! 


(1-15) 


This  relation  defines  the  diagonal  elements  of  CM 


<CM> 


ii 


>ii 


(1-16) 


The  off-diagonal  elements  of  Cy  are  determined  by  the 
following  prescription 


Qj  -  -Qi  =  Q 

Qk  =  0,  k  f  i,j  (1-17) 
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With  this  assumption,  we  can  expand  Equation  ( I — i 4 )  for 
i  and  j 


vj  *  9  -  9 

v±  =  (C-1)  Q  -  <c_1)ii  Q  (1-18) 

So  that 

V.-V.  =  /(C_1)  .  .  +  (C~1) .  . -(C-1) . .-(C-1) .  . )  Q  (1-19) 
J  1  'jj  11  Ji  ij  / 

Re-writing,  we  have 

Q  =  (CM)ij  AV 


(CM)  =  — -I - - — - 

J  (C  ) ■  ■  +  (C  x)  .  .  -  2 ( C  x) .  . 
^  '11  v  '.JJ  v 


AV  =  Vj  -  Vi 


(1-20) 


The  manipulations  outlined  above  result  in  a 
capacitance  matrix  C^j,  the  elements  of  which  can  be 
conveniently  measured.  In  this  definition  of  Cy,  (C^)^ 
is  the  ratio  of  the  charge  on  the  i'th  conductor  to  the 
potential  on  the  i'th  conductor  when  all  other  conduc¬ 
tors  have  no  charge.  (CM). .  (i^j)  is  the  ratio  of  the 

vJ 

charge  on  the  i'th  conductor  to  the  potential  difference 
Vj-Vi  when  all  other  conductors  k(k^i,j)  have  no  charge. 


Appendix  II 

SOLUTION  OF  TRANSMISSION-LINE  EQUATIONS 


In  this  appendix  we  will  discuss  the  solution 
to  the  T-L  equations  as  presented  in  the  main  body  of 
the  text.  In  order  to  arrive  at  the  exact  solution  we 
have  made  the  following  assumptions: 


1)  There  is  only  one  propagation  velocity,  V. 

2)  Iq  is  uniform  along  the  multiwire  bundle 
(normal  X-ray  incidence). 

3)  Resistive  terminations. 


The  equations  to  be  solved  are: 


3V  _  =  31 

3X  3t 


31 

ax 


+  i 


D 


(II-l) 


together  with  the  boundary  conditions 


I(X,t=0)  =  V(X,t=0)  =  0 


V(X=0,t)  =  -ZL  I(X=0,t) 


V(X=  5. ,  t )  =  Zn  I(X-t.t) 


( 1 1  -  2) 
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where  Zl(Z^)  is  the  termination  impedance  matrix  on  the 
left  (right)  end  of  the  cable  of  length 


The  method  of  solution  is  to  take  the  fourier 
transform  of  the  two  equations  using  the  following  con¬ 
vention 


i(X,  w) 


V(X, w) 


, .  „  .  .  iwt 

dt  I  ( X ,  t )  e 


dt  V(X , t )elwt 


f(w) 

D 


f  dt  ID(t)eiwt 

J0 


(  1 1 -3  ) 


The  fourier  transformed  version  of  these  two  equations 
are 


3 V 
3X 


iwL  I 


3? 

3X 


iwC  V  +  ID 


(II-4) 


We  can  take  the  second  derivative  of  each  of  these 
equations  with  respect  to  X,  noting  by  assumption  1 
that  CL=LC=^  T  and  by  assumption  2  that  1^.  =  0, 

and  arrive  at  the  following  two  equations 
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1L 


•  * 


0 


32I 

ax2 


2— 
+  k^I 


2- 

afv 

ax2 


2— 
+  K  V 


=  iwL  I, 


where 


2  w 


The  solution  of  these  two  equations 


is 


~  =  eikX  y  +  e'ikX  v  + 

4" 


iZ0ID 


f  ,  eiM  T+  ♦  e-ikX  T 


where 


f  _  KjcJb 

zo  v 

v+  =  z0  7+ 


V_  =  Z0  X- 


=  =  =  2ik  J,. -1 
I+=<1-PRPL6  } 


<W’1'5 


iZ0ID 
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(II-5) 

( II-6 ) 

straightforward 

(II-7) 


ir». 


—  =  =  =  2ik  1  ^  =  =  -i  =  ====  _i 

I-=  ‘^'Vl6  >  <VZR>  -  VW  e  Jx 


Rv“0  L' 


ikA  I  lZ0rD 


PR  ‘  <z0  +  V1  <z0  -  Sr) 


PL  *  <Z0  *  ZLrl  <Z0  -  ZL> 


(11-8) 


The  solution  for  I  at  X=0  is  simply 


I(X-0,w)  = 


K  +  X- 


*  (<"-PL>[1-PRPLe2ikl 


-if  -  2iki 

VZ0+V  e 


+  (ZQ+2R)~1elk£  -  (Z^Z^*1!  ( 1 1  -9 ) 


We  can  take  the  inverse  fourier  transform  of  this  equa¬ 
tion,  according  to  the  following  definition 


-  2 kf. 


dw  I  ( X ,  w )  e 


V(X,t) 


2ttJ 


dw  V(X,w)e 


ID(t) 


-  -Lf 

2n  J 


dw  I^(w)e 


-iwt 


(11-10) 
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The  analysis  of  the  fourier-transform  of  Equation  (II-9) 
is  made  intelligible  by  the  following  observations.  Let 
Q(t)  be  the  total  charge  deposited  at  time  t  due  to  the 
Norton  equivalent  driver  1^ 


-f* 


Q(t)  =  J  In(Odt 
0 


(11-11) 


Using  Equation  (II- 10)  we  can  write 


Q(t 


/t  /-+  °° 

dt'  J  dw  fD( 


w)e 


-iwt 


„_i  ( 
2-n  J 


+°°  ilD(w)e 

dw  - 

w 


-iwt 


(11-12) 


So  that  we  have 


if 

2 7T  I 


+oo 


dw  e 


-iwt 


=  V  ZQ  Q(t) 


(H-13) 


In  addition  we  have 


-  -  2iki  _1  »  =  -  n  i2nk£ 

^-PRPLe  >  "  Z  <PRPL>  e 

n=0 


(11-14) 
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I 


So  that 


1(0, w)  = 


(1-PL)  f  (PRPL)Qe 
ln=0 


i2nkJ l 


2ik£  —  —  ikS. 

P^I^^)'1®  +  CZQ+ZR)  e 

«  =  -4  1  vd] 

-  cvV  )  HH 


If  we  take  the  fourier-transform  of  this  equation, 
using  Equation  (11-10) ,  we  have 

1(0, t)  =  -  |  (T  +  PL)  x 


(11-15) 


£  (PRPL)n  |«Q(t-2nT£)  -  PR  6Q(t-(2n+l)t£)|  (II  16) 


where 


5Q( t )  =  Q(t)  -  QCt-T^) 


-  K  / 

Tl  ~  J V 


(11-17) 


as  the  final  solution  for  the  current  that  flows  at  X-0  as 


a  function  of  time. 


Appendix  III 


ENERGY  THEOREM 


An  important  consequence  of  the  properties  of 
the  T-L  equations  as  presented  in  the  main  body  of  the 
text  and  in  the  preceding  Appendix  is  the  presence  of 
an  energy  theorem.  The  total  energy  density  for  a 
multiwire  bundle  is  defined  by  the  following  equation 


U 


i  r.  v 


( V , C  V  )  +  (I,  LI) 


(III-l) 


where  the  inner  product  of  two  vectors  A  and  B  is  denoted 
by  the  symbol  (A,B).  If  we  take  the  partial  derivative 
of  y  with  respect  to  t  and  X  and  use  the  T-L  equations, 
Equation  ( I I — 1 )  to  simplify,  we  can  write  a  Poynting 
equation  for  y 


ajj  +  as  = 

3t  3X 


,V) 


(III-  2) 


where  S  is  the  energy  flux  (Poynting  Vector)  and  is 
defined 


S  =  (I,V) 


(  1 1 1  -  3) 


and  the  source 
transported  to 
term,  which  is 
drivers  Tq. 


term  is  (Ip.V).  The  energy  which  is 
the  loads  is  generated  by  the  source 
directly  related  to  the  Norton  Equivalent 
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If  we  integrate  Equation  ( 1 1  —  2 )  over  the 
length  of  the  cable  l,  we  obtain  an  expression  for 
the  total  energy  in  the  multiwire  cable  E 

+  |(KO,  V(0))  -  (1(0),  7(0)) 

l 

(ID,  V)  dx  ( I II- 4) 


If  we  again  assume,  as  we  did  in  Appendix  II,  that  Ip  is 
independent  of  X,  then  the  right-hand  side  of  the 
previous  equation  can  be  written 

l 

dx  (ID,  V)  =  £(Id,  <  V »  (HI-5) 


where  <V>  is  the  average  voltage 


<V> 


£ 

V  dx 


The  quantity,'  therefore,  in  Equation  ( 1 1 1  —  5 )  represents 
the  total  power  source  and  the  total  available  energy  is 


e  -  if  dD,  <v»  dt 

Jo 


( II 1-6) 


We  can,  again,  use  the  T-L  equations  to  get  an  expression 
for 
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15T<V>  -  C’1  ID  -  (l(t)  -  1(0))  UII-7) 

This  equation,  together  with  the  relationship  between 
current  driver  and  charge 


allows  us  to  rewrite  the  total  available  energy  as 


1(0) 


>i 


1 

2 


(  1 1 1  -9  ) 


where  QD  is  the  total  charge  transfer  in  the  multiwire 
cable 


QD 


(HI  -10  ) 


To  arrive  at  the  expression  in  Equation  (III-9)  for  S, 
we  have  assumed  that  at  t  =  <*>,  V(<*>)=(5. 

Since  the  expression  for  £  in  Equation  (III-9) 
represents  the  total  energy  that  can  be  dissipated  across 
all  terminations  in  a  multiwire  bundle,  it  is  of  interest 
to  cast  this  equation  into  a  more  useable  form.  It  is 
immediately  obvious  that  the  energy  is  not  a  constant  of 
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the  motion,  since  it  depends  on  the  solution  1(1),  1(0), 
which  is  associated  with  the  fact  that  the  Poynting 
equation,  Equat ion  (  1 1  I  -  3 )  has  a  current  source  term. 

We  can,  however,  arrive  at  a  usable  expression  for  £, 
since  we  have  an  analytic  solution  for  the  current  T 
that  flows  across  the  terminations,  as  presented  in 
Appendix  II.  For  convenience  in  this  derivation,  we 
will  assume  that  the  terminations  are  symmetric,  i.e. 

5L  =  ?R  =  F  (HI-11) 

The  solution  for  the  current  at  X=0  is,  therefore 

1(0, t)  =  -  |  (f  +  F) 


£(P)2n  [6Q(t-2nTA)  -  P  6Q(t-(2n+l)xJl)]  (III-12) 

For  the  case  of  symmetric  terminations,  we  have 

1(0  =  -  7(0)  (ill -13 ) 

So  that 
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In  addition,  we  will  assume  that  the  waveform  f(t)  for 
the  current  driver  7^  is  a  square  wave  of  width  Xp,  so 
that 


IdCO 

=  l0  f(t) 

Q(t) 

=  v  0  1*  £  Tp 

and 

■  Totp  l±TP 

qd 

T0tp 

(III -15 ) 

where  Iq  is 

the  peak  driver  in  a/m. 

The  other  time  of 

interest  is 

the  propagation  time,  x^ 

where 

II 

£/ V 

(III-16) 

In  order  to  obtain  an  expression  for  £  from  Equation  (III-14) 
we  must  assume  a  particular  value  for  the  ratio 


^  =  Tp/T£ 


(III-17) 


Since  we  are  interested  in  obtaining  an  expression  for  £ 
the  maximum  energy  that  can  be  dissipated,  it  is  natural 
to  look  at  the  particular  values  of  £  that  correspond  to 
long  cables,  i.e. 


0  <  c  <  1 


(III -18 ) 


If  we  assume  the  range  of  <;  as  given  above,  and  substi¬ 
tute  the  expression  for  1(0)  into  Equation  (III-14), 
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y 


« 


using  Equation  (III- 15)  we  can  perform  the  integrations 
over  time  and  sum  the  infinite  series.  The  result  is  a 
closed  form  expression  for  the  total  energy  that  can  be 
dissipated  across  all  terminations 


£/<?max  -  1  -  |  Cd  +  <P>)  (III- 19) 

where 


This  result  is  valid  in  the  range 


0  <  c  <  1  ( I 11-21 ) 

In  general,  the  solution  for  S  is  different  in  other 
ranges  of  5,  getting  more  complicated  as  £  gets  larger. 

The  result  for  the  range 


1  <  5  <  2 


can  be  written 


( III-22) 


<?/<?MAX  =  1  “  [t  (1+2  <P>  +  <P2>) 

-  |  C  (<P>  +  <P2>) J  1  £  C  £  2 


(HI-23) 


t 
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where 


Appendix  IV 


1 


1  ! 


MULTIWIRE  TERMINATION  IMPEDANCE 
MATRIX  AND  REFLECTION  COEFFICIENT 

The  solutions  to  the  T-L  equations  presented 
previously,  both  for  the  current  and  energy  into  the 
loads,  depend  on  the  terminations  through  the  impedance 
matrix  Z  or,  alternatively,  on  the  reflection  coeffi¬ 
cient  p  which  is  a  function  of  Z.  In  this  appendix 
we  will  consider  the  relationship  between  Z  and  physical 
resistance  and,  also,  what  the  reflection  coefficient  P 
looks  like  for  the  various  possible  terminations. 

IV. 1  IMPEDANCE  MATRIX 

The  boundary  conditions  for  an  n+1  bundle  of 
wires  terminated  at  the  left  (right)  end  with  a  termina¬ 
tion  impedance  matrix  Z l  (Zr)  is: 

V(0,t)  -  -  ZL  1(0, t) 

V(i,t)  =  fR  I(£,t)  (IV-1) 


If  we  consider  the  left  side  of  the  transmission  line, 
we  can  write: 


I 


9 


If  we  denote  the  current  (voltage)  at  X=0  in  wire  i  by 
1^,  (V,)  we  have 


Xi 


N 

-E 

j=i 


(Z71) . .  v . 

L  'ij  J 


(IV-3) 


This  current  will  divide  depending  upon  the  physical 
resistances  that  it  sees.  If  we  denote  these  physical 
resistances  by  Rjj  j=l,  n,  where  Rii  is  the  resistance 
from  i  to  ground  and  Rjj  ,  j=l,  n  ( i^ j )  are  the  resis¬ 
tances  from  i  to  all  the  other  wires  in  the  bundle, 
then  we  have 


( IV-4  ) 


Equating  these  two  expressions  for  Ii(  we  have 


N  vi  -  V  V  N  , 

E  - 1  +  =  E  <ZT,  >n  V. 

j=l 


Rij 


Rii  A  '“L  j 


(IV-5) 


Rearranging  these  terms  we  have 


N 

j=l  Rij 


*  sil  -  <2")- 


V  . 

l 


N 

-  E 

j=i 


£  * 


V.  =  0 


(IV-€  ) 


» 
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These  equations  imply 


(Z71) .  . 
v  L  n 


N  i 

T  — 

R.  . 

J  =  1 
if  3 


R-  Rii 


(z:1) . . 

L  'ij 


Rij 


We  can  rewrite  the  first  equation  using  the  second 
equation  and  we  obtain 


(ZT 1) . . 
L  -'ll 


N 

E 

j=i 

ifj 


(Z.1)  .  .  +  ~ 

L  'ij  R_.  j 


11 


which  implies 


Rii 


<ZLhi  *  &  ‘Oij 

ifj 


N 


.-1, 


■  £  (ZL  >ij 


And  so  we  obtain  the  following  two  equations 


aij 


i^j 


Rii 


E  (ZL1}ij 

j  =  l 


( IV-7  ) 


( IV-8  ) 


(iv-9 ; 


(IV-10) 
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If  we  desire  to  characteristically  terminate  the  multi¬ 
wire  cable,  we  can  use  the  above  two  equations  to 
determine  the  appropriate  physical  resistances.  If  we 
desire  to  determine  what  termination  matrix  Z l  to  use 
corresponding  to  a  given  set  of  physical  resistances 
Ri j ,  we  can  use  Equation  (IV-7).  These  can  be  written 
as  : 


>ii 


(IV-11) 


IV. 2  EXAMPLES  -  TWO  WIRE 


In  this  section  we  will  consider  some  particu¬ 
lar  examples  of  the  results  indicated  in  the  previous 
paragraph.  Let  us  consider  a  TSP  with  the  following 
capacitance  C,  inductance  L  and  characteristic  impedance 
matrix  ZQ : 


207. 

-  36. 

36. 

207. 

223. 

11. 

11. 

223. 

33.3 

3. 

3.73 

33. 

pf/ra 


nH/m 


( IV- 12) 
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CASE  1  : 


CHARACTERISTIC  TERMINATION 


If  we  desire  to  terminate  the  TSP  character¬ 
istically,  we  can  use  Equation  (IV-10)  to  determine  the 
appropriate  physical  resistances.  The  inverse  of  the 
characteristic  impedance  matrix  is  given  by: 


3.04 

x  10  2 

-  3.41  x 

io-3\ 

> 

1 

3.41 

x  10-3 

3.04  x 

10-2/ 

(IV 

and  the  physical  resistance  matrix  R  is: 


(32.9 
293.6 


293.6 

32.9 


A 


( IV-14 ) 


The  reflection  coefficient  p  when  the  cable  is  charac¬ 
teristically  terminated  is  given  by: 


P 


0 


o) 


(IV-15) 


CASE  2:  WIRE-TO-SHI ELD  TERMINATION 

If  we  desire  to  terminate  the  TSP  only  wire-to- 
shield  with  Ri2=^21=00>  corresponding  to  no  wire-to-wire 
termination,  we  can  set  Ri2=^21  =  0°  •'•n  Equation  (IV-11)  to 
determine  the  appropriate  Z,  matrix: 
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* 


<ZZ\l  -  E 


11 


(z;1) 


L  '22  R 


22 


*  «Z£1>21  ■  0 


(IV- 16) 


which  implies: 


R 


11 


22 


(IV-17  ) 


For  this  case,  therefore,  the  diagonal  elements  of  Z^ 
are  identical  to  the  physical  resistances  ,  R22‘ 

If  we  write  the  characteristic  impedance 
matrix  in  the  following  way 


Z 

Z, 


( IV-18  ) 


then  p  can  be  written: 


P  “  K 


(Z1-R1)(Z2+R2)-Zz  2ZR2 


2ZR, 


(Z1+R1)(Z2-R2)-Z< 


D  =  (Z1+R1)(Z2+R2)-Zi 


( IV-19  ) 
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9 


% 


It  is  important  to  note  from  the  above  expression  for 
p  that  whenever  the  two  wires  are  terminated  wire-to- 
shield  only,  p  can  never  be  diagonal  for  any  finite 
value  of  R^  and  R2* 


CASE  3:  SHORT-CIRCUIT  WIRE-TO-WIRE 


If  we  desire  to  determine  the  appropriate  Z 
when  the  two  wires  of  the  TSP  are  tied  together  and 
terminated  with  one  resistance  to  ground  we  can  set: 


12 


=  R 


21 


=  0 


lll 


=  R 


22 


=  R 


This  corresponds  to: 


(IV-20) 


(IV-  21) 


The  reflection  coefficient  is: 


P 


1 

D 


/ (Z1-R)(Z2+R)-Z2 
\2ZR 


2ZR 

(Z1+R)(Z2-R)-Z2 


D  =  (Z1+R)(Z2+R)-Z2 


(IV- 22) 
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For  the  general  situation  of  any  combination  of 
physical  resistances  we  can  write  down  the  appropriate 
termination  matrix  Z^. 


R  = 


ZL  = 


'll 

R12 

'21 

R22 

1 

/ 

(R22+R21)(R11R12)  R12R11R22 


12  21  11  21  22  12\r2iRi1R22  (R^+R^)  <*21R22>  J  (IV‘23) 


Appendix  V 


TWO  WIRE  SOLUTIONS 


V.1  CURRENT 

The  termination  dependence  of  both  the  current 
and  energy  in  a  multiwire  bundle  can  be  conveniently 
represented  by  the  dimensionless  matrix  p  defined  as 

P  -  <ZQ  +  Z)'1  (!0  -  f)  ( V-l ) 


The  exact  solution  for  the  current  given  in  Equation 
(11-16)  indicates  that  the  fact  that  p  is,  in  general,  a 
non-diagonal  matrix  distinguishes  a  multiwire  bundle  from 
a  single  wire  transmission  line.  Since  the  characteristic 
impedance  matrix  Zq  can  never  be  diagonal  for  any  physi¬ 
cally  reasonable  multiwire  bundle,  p  can  never  be  diagonal 
for  any  finite  value  of  the  terminations.  It  is  possible, 
however,  to  make  a  transformation  which  diagonalizes  p,  and 
which  allows  the  solution  for  a  multiwire  bundle  to  be 
written  in  terms  of  single  wire  solutions.  While  this 
procedure,  in  principle,  can  be  done  for  any  number  of 
wires,  in  practice,  the  algebra  becomes  rather  tedious 
for  n>2.  In  the  following  paragraph  we  will  sketch  this 
procedure  for  a  two-wire  bundle. 


The  characteristic  impedance  matrix  is  defined 


(V-2  ) 
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Z 


and  the  termination  impedance  matrix  is  defined 


( V-3  ) 


where  we  will  only  consider  terminations  wire-to-shield, 
and  R  and  R  are  the  physical  resistances  wire-to- 
shield.  The  reflection  coefficient  can  be  written 


P  =  A 


(Z1-R1)(Z2+R2)-Z^  2ZR2 


2ZR1  (Z1+R1)(Z2-R2)-Z^y 


where 


A  =  (Z1+R1)(Z2+R2)-Z" 


( V-4) 


The  e-values  of  p  can  be  found  directly 

2  (l+r1)(l+r2)-C2 

where 

rl  *  Rl/Zl 
r2  ~  R2^22 
C  =  Z/CZ^)1/2 


(V-5) 


(V-6) 


IPO 


The  e-vectors  corresponding  to  these  e-values  can  also 
be  obtained  directly 


(8+ct)/(l+OB+(l-s)a 

?(8+a)/(l+06+(l-i;  >a 


1/2 


1/2/ 


(V-  7  ) 


where 


a  -  R^Z2  -  ^2^1 
6  =  (a2  +  4Z2  R1R2)1/2 

C  =  R1/R2  (V- 8  ) 

The  quantity  C  is  a  measure  of  the  wire-to- 
wire  coupling.  It  is  easily  seen  that  the  solution 
reduces  to  two  independent  uncoupled  transmission 
lines  when  C  is  zero.  From  general  properties  of  the 
impedance  matrix,  it  may  be  shown  that  the  maximum 
value  of  C  is  unity. 

Once  the  e-values  and  e-vectors  have  been 
determined,  it  is  straightforward  to  construct  the 
transformation  matrix  T.  If  we  denote  e^  and  e2  as 
the  unit  vectors  in  the  original  basis,  which  can  be 
written  as 
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I 


We  can  write 


In  order  to  write  the  solution  for  the  current  in  wire 
one  and  two  it  is  necessary  to  obtain  (T)"1  which  we 
can  write: 


li<2 


and 


D  =  (l+oV  -  (l-02a2 


(V-12) 


In  order  to  understand  wire-to-wire  coupling 
in  a  multiwire  transmission  line,  it  is  of  interest  to 
drive  one  wire  and  look  at  the  current  that  flows  in 
the  other  wire.  The  analysis  that  we  have  just  developed 
allows  us  to  accomplish  this.  Let  X  denote  the  basic 
current  driver 

X  =  XiSi  (V-13) 


From  the  previous  section  we  have 


(T  1)  .  . 


n  . 
J 


So  that 


X  =  (T-1)ij  a.  (V-i4) 

The  current  that  flows  in  both  lines  can  then  be  written 

h  -  I(V  V  (V-15) 
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where  I(Aj)  is  the  single  wire  solution  for  the  e-value 
A j .  For  the  special  case  that  we  drive  only  wire  two, 
we  have 


xx  =  0 


x2  =  1 


(V-16) 


and  the  solution  for  the  current  in  wire  one  and  two  can 
be  written: 


/  2  2  \  1/2 
■■  •  te?) 


I  (  A,  )  -  I(A„) 


r2  '  (^f)  +  (‘if)  ‘‘V 


(V-17) 


where  I(A2)  are  the  single  wire  solutions  with  reflec¬ 


tion  coefficient  X2- 


V.2  ENERGY 

Once  the  current  in  both  lines  is  written  in 
terms  of  single  wire  solutions,  it  is  of  interest  to 
compute  the  energy  dissipated  across  the  terminations. 
Our  main  interest  in  this  analysis  is  the  calculation 
of  the  energy  dissipated  across  termination  one  when 
wire  two  is  driven,  since  all  of  this  energy  represents 
coupled  energy.  The  current  that  flows  in  wire  one  can 
be  written 


ZR, 


I-.  = 


[(RlZ2-R2Zl)2t4z2RlE2] 


1/2 


jl(X1)-I(X2)J  ( V-18) 


and  the  energy  is 


E1  =  R1A! 


CO 

U  dt  [l(X1)2+I(X2)2-2I(X1)I(X2) 


where 


ZR, 


A  = 


[(R1Z2-R2Z1)2+4Z2R1R2 


il/2 


( V-19 ) 


If  we  assume  a  square  wave  driver,  we  can  perform  these 
integrations  as  has  been  done  in  previous  sections.  In 
particular,  if  we  take  the  long  line  limit  which  corres¬ 
ponds  to 


C  =  VT£  <<<A 


( V-20 ) 


We  arrive  at  the  following  expression  for  E ^ : 


E1  =  E2 


MAX 


1 

2 


r  r  2C2 
rlr2  ^ 


(r'i-r2)2+4c2rir2 


1+xi  i+x2 
1~X  i  l-*2 


-  2 


(1+X1)(I>X2) 

(1-A1x2) 


(  V-21  ) 


» 
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MAX 

Where  E2  is  the  maximum  energy  capability  of  wire  2 
and  is  given  by 

e2MAX  =  |(VZ2)QD2£  ( V-22 ) 

From  these  expressions,  one  can  show  that  the  coupled 
current  and  energy  tend  to  zero  as  the  termination 
resistance  on  line  2  goes  to  zero  (shorted)  and 
increases  as  the  termination  resistance  is  increased. 

If  line  2  is  opened,  the  current  on  line  1  simply 
appears  as  if  line  1  were  being  directly  driven  with 
an  effective  source  term  of  C  times  the  line  2  source 
term.  In  particular,  all  the  energy  capability  of 
line  2  is  not  coupled  into  line  1.  The  final  voltage 
on  the  second  line  is  not  zero,  and  leads  to  capaci¬ 
tive  energy  storage  of  a  portion  of  the  energy  capa¬ 
bility  of  this  line. 
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Appendix  VI 


TWO-WIRE  CALCULATIONS 


We  have  considered  a  number  of  two-wire  cables. 
Each  cable  is  labeled  by  the  designation  TWI  where 
1=1,2, 3.  The  following  parameters  distinguish  each 
cable  from  one  another. 


TWI 


207. 


-36.0 


TW2 


7, 


36.0 

207 

223. 

11 

11. 

223 

33 . 3 

3 

3.73 

33 

207. 

-  36 

36. 

207 

55.4 

9 

pf/m 


nH/m 


ohms 


pf/m 


nh/m 


ohms 


TW3 


226. 

U109. 


-109. 

751. 


pf/ra 


L 


8. 

16. 


nH/ra 


(15.9 
23.5 


23.5 

4.8 


ohms 


In  the  following  pages  we  will  present  the  results  of  the 
calculations.  The  cable  designated  TW1  is  the  two-wire 
cable  which  has  two  phase  velocities  whereas  cables  1  and 
3  each  have  one  phase  velocity.  Each  plot  presented  is 
labeled  by  a  header  which  indicates  the  parameters  used. 
In  the  following  Tables  we  indicate  the  definition  of 
these  headers.  For  cables  TW1  and  TW2  the  current 
drivers  were  as  follows 


and  for  TW3  the  current  drivers  were 


1 


TABLE  VI. 1 


TW1 


33.3 

0 


R  =  termination  impedance  matrix  (ohms)  (symmetric 
terminat ion ) 

Tp  =  pulse  width  (square  wave)  (ns) 

L  =  length  of  cable  (meter) 


TP 

L 

Header 

1Z 

3 

.  3 

TW1R11T1L1 

1Z 

3 

3. 

TW1R11T1L2 

1Z 

3 

30. 

TW1R11T1L3 

1Z 

30 

.3 

TW1R11T2L1 

1Z 

30 

3. 

TW1R11T2L2 

1Z 

30 

30. 

TW1R11T2L3 

1Z 

300 

.3 

TW1R11T3L1 

1Z 

300 

3. 

TW1R11T3L1 

1Z 

300 

30. 

TW1R11T3L3 

z 

3 

.3 

TW1R22T1L1 

z 

3 

3. 

TW1R22T1L2 

z 

3 

30. 

TW1R22T1L3 

z 

30 

.  3 

TW1R22T2L1 

z 

30 

3. 

TW1R22T2L2 

z 

30 

30. 

TW1R22T2L3 

0 

33.3, 


ohms 
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TABLE  VI. 1  (Continued) 


R 

TP 

L 

Header 

Z 

300 

.3 

TW1R22T3L1 

2 

300 

3. 

TW1R22T3L2 

2 

300 

30. 

TW1R22T3L3 

10  2 

3 

.3 

TW1R33T1L1 

10  2 

3 

3. 

TW1R33T1L2 

10  2 

3 

30. 

TW1R33T1L3 

10  2 

30 

.3 

TW1R33T2L1 

10  2 

30 

3. 

TW1R33T2L2 

10  Z 

30 

30. 

TW1R33T2L3 

10  Z 

300 

.3 

TW1R33T3L1 

10  Z 

300 

3. 

TW1R33T3L2 

10  Z 

300 

30. 

TW1R33T3L3 

TABLE  VI. 2 


TW2 


R  =  termination  impedance  matrix  (ohms)  (symmetric 
termination) 

Tp  =  pulse  width  (square  wave)  (ns) 

L  =  length  of  cable  (meter) 


R 

TP 

L 

Header 

.  1Z 

3 

.3 

TW2R11T1L1 

•  1Z 

3 

3. 

TW2R11T1L2 

.  1Z 

3 

30. 

TW2R11T1L3 

.1Z 

30 

.3 

TW2R11T2L1 

•  1Z 

30 

3. 

TW2R11T2L2 

.  1Z 

30 

30. 

TW2R11T2L3 

•  1Z 

300 

.3 

TW2R11T3L1 

.1Z 

300 

3. 

TW2R11T3L2 

.  1Z 

300 

30. 

TW2R11T3L3 

z 

3 

.3 

TW2R22T1L1 

z 

3 

3. 

TW2R22T1L2 

z 

3 

30. 

TW2R22T1L3 

z 

30 

.  3 

TW2R22T2L1 

z 

30 

3. 

TW2R22T2L2 

z 

30 

30. 

TW2R22T2L3 
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TABLE  VI. 2  (Continued) 


R 

TP 

L 

Header 

z 

300 

.3 

TW2R22T3L1 

z 

300 

3. 

TW2R22T3L2 

z 

300 

30. 

TW2R22T3L3 

10 

z 

3 

.3 

TW2R33T1L1 

10 

z 

3 

3  . 

TW2R33T1L2 

10 

z 

3 

30. 

TW2R33T1L3 

10 

z 

30 

.3 

TW2R33T2L1 

10 

z 

30 

3. 

TW2R33T2L2 

10 

z 

30 

30. 

TW2R33T2L3 

10 

z 

300 

.3 

TW2R33T3L1 

10 

z 

300 

3. 

TW2R33T3L2 

10 

z 

300 

30. 

TW2R33T3L3 

.1Z 

1+Z2 

3 

3 

TW2R12T1L2 

.  1Z 

1+10Z2 

3 

3 

TW2R13T1L2 

V 

•1Z2 

3 

3 

TV72R2 1T1L2 

Z1+10Z2 

3 

3 

TW2R23T1L2 

10  z 

1+-1Z2 

3 

3 

TW2R31T1L2 

10  z 

1+Z2 

3 

3 

TW2R32T1L2 
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TABLE  VI. 3 


TW3 


R  =  termination  impedance  matrix  (ohms)  (symmetric 
termination) 

Tp  -  pulse  width  (square  wave)  (ns) 

L  =  length  of  cable  (meters) 


R 

T 

p 

L 

He ader 

.  1Z 

3 

3 

TW3R11T1L2 

Z 

3 

3 

TW3R22T1L2 

10  Z 

3 

3 

TV/3R33T1L2 

•  1Z1+Z2 

3 

3 

TW3R12T1L2 

.  lZ^lOZg 

3 

3 

TW3R13T1L2 

Zl+ • 1Z2 

3 

3 

TW3R21T1L2 

z1+ioz2 

3 

3 

TW3R23T1L2 

ioz^.  iz2 

3 

3 

TW3R31T1L2 

lOZ^Zg 

3 

3 

TW3R32T1L2 
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Appendix  VII 


SEVEN-WIRE  CALCULATIONS 


We  have  considered  various  distributions  of 
drivers  for  the  seven  wire  bundle  discussed  in  Section  V. 
These  distributions  represent  typical  groups  of  seven 
wires  in  the  thirty-seven  wire  bundle  also  discussed  in 
Section  V.  The  propagation  time  and  pulse  width  were 
chosen  such  that 


which  corresponds  to  a  length  of  cable  equal  to  ^3m. 
Three  different  distributions  of  drivers  were  considered 
SE1 ,  SE2 ,  SE3  according  to  the  following 
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Each  plot  is  labeled  by  a  header.  The  definition  of  the 
header  for  each  plot  is  given  in  Table  (VII . 1)-(VII . 3 ) . 
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TABLE  VII. 1  (Continued) 


R  =  termination  impedance  matrix  (ohms)  (symmetric 
termination ) 

Tp  =»  pulse  width  (square  wave)  (ns) 

L  =  length  of  cable  (meter) 


R 

Tp 

L 

Header 

.  1Z 

17 

3 

SE1R11T1L1 

Z 

17 

3 

SE1R22T1L1 

10  Z 

17 

3 

SE1R33T1L1 

. lZj+Zg 

17 

3 

SE1R12T1L1 

.  IZj+IOZjj 

17 

3 

SE1R13T1L1 

Z1+’1Z2 

17 

3 

SE1R21T1L1 

z1+ioz2 

17 

3 

SE1R2 3T1L1 

lOZ^.  1Z2 

17 

3 

SE1R31T1L1 

lOZj+Zg 

17 

3 

SE1R32T1L1 
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TABLE  VII. 2  (Continued) 


R  =  termination  impedance  matrix  (ohms)  (symmetric 
termination) 

tp  =  pulse  width  (square  wave)  (ns) 

L  =  length  of  cable  (meter) 


R 

TP 

L 

Header 

.  1Z 

17 

3 

SE2R11T1L1 

Z 

17 

3 

SE2R22T1L1 

10  z 

17 

3 

SE2R33T1L1 

•lzi+z2 

17 

3 

SE2R12T1L1 

.1Z1+10Z2 

17 

3 

SE2R13T1L1 

Z1+'1Z2 

17 

3 

SE2R21T1L1 

V10Z2 

17 

3 

SE2R23T1L1 

10Z1+.1Z2 

17 

3 

SE2R31T1L1 

1°Z1+Z2 

17 

3 

SE2R32T1L1 
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TABLE  VII. 3  (Continued) 


R  =  termination  impedance  matrix  (ohms)  (symmetric 
termination ) 

Tp  =  pulse  width  (square  wave)  (ns) 

L  =  length  of  cable  (meter) 


R 

TP 

L 

Header 

.  1Z 

17 

3 

SE3R11T1L1 

z 

17 

3 

SE3R22T1L1 

10  z 

17 

3 

SE3R33T1L1 

•1Z1+Z2 

17 

3 

SE3R12T1L1 

.1Z1+10Z2 

17 

3 

SE3R13T1L1 

Zl+ . 1Z2 

17 

3 

SE3R21T1L1 

z1+ioz2 

17 

3 

SE3R22T1L1 

ioz1+.iz2 

17 

3 

SE3R31T1L1 

10Z1+Z2 

17 

3 

SE3R32T1L1 
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Appendix  VIII 


LOSSY  CABLE  CALCULATIONS 


The  effects  of  loss  have  been  considered  for 
the  following  cable:' 

RG  -  316 
C  =  95  pf/m 

L  =  237.5  nH/m 

e  =  2.05 

In  the  following  pages  we  will  list  the  results  of  the 
calculations.  Each  plot  is  labeled  by  a  header  to 
identify  it.  In  Table  (VIII. 1)  is  indicated  the  defi¬ 
nition  of  the  header.  In  the  first  two  plots  are 
indicated  R(W)  and  G(W)  as  is  used  in  all  the  calcula¬ 
tions  and  the  remainder  of  the  plots  indicate  the 
results  of  the  calculations. 
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TABLE  VIII. 1 


R  =  termination  impedance  (ohms)  (symmetric 
termination) 

Tp  =  pulse  width  (square  wave)  (ns) 

L  =  length  of  cable  (meters) 


R 

TP 

L 

Header 

.5 

3 

.3 

0N1R1T1L11 

.5 

3 

3. 

ON1R1T1L21 

.5 

3 

30. 

ON1R1T1L31 

.5 

30 

.3 

ON1R1T2L11 

.5 

30 

3. 

ON1R1T2L21 

.5 

30 

30. 

ON1R1T2L31 

.5 

300 

.3 

ON1R1T3L11 

.5 

300 

3. 

ON1R1T3L21 

.5 

300 

30. 

ON1R1T3L31 

50 

3 

.3 

ON1R2T1L11 

50 

3 

3. 

ON1R2T1L21 

50 

3 

30. 

0N1R2T1L31 

50 

30 

.  3 

ON1R2T2L11 

50 

30 

3. 

0N1R2T2L21 

50 

30 

30. 

ON1R2T2L31 

50 

300 

.  3 

ON1R2T3L11 

50 

300 

3. 

0N1R2T3L21 

50 

300 

30. 

0N1R2T3L31 

5000 

3 

.3 

0N1R3T1L11 

5000 

3 

3. 

ON1R3T1L21 

5000 

3 

30. 

ON1R3T1L31 
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TABLE  VIII. 1  (Continued) 


R 

TP 

L 

Header 

5000 

30 

.3 

ON1R3T2L11 

5000 

30 

3. 

0N1R3T2L21 

5000 

30 

30. 

ON1R3T2L31 

5000 

300 

.3 

ON1R3T3L11 

5000 

300 

3. 

ON1R3T3L21 

5000 

300 

30. 

ON1R3T3L31 
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LIST  OF  TERMS 


SGEMP 

DSCS  III 

C.I.T. 

L/C 

TL 

2 

2c 

2x 

i'o(x>+) 

r 

* 

v(x) 

0 

Q 

*D 

S,1 

2 

Il.(l2) 


IL 

TP 

voc 


System  Generated  Electromagnetic  Pulse 

Defense  Satellite  Communications  System,  Phase  III 

Current  Injection  Test 

Inductive/Capacitance 

Transmission  Line 

Impedance 

Characteristic  line  impedance 

Load  impedance  at  position  x 

point  current  source 

reflection  coefficient 

line  propagation  constant 

line  voltage  at  position  x 

angle  of  incidence  of  photon  to  cable 

total  charge  accreted  on  open  circuited  line 

Current  Driver 

line  length 

distance  of  line  from  its  end 

contribution  to  current  from  all  source  left 
(right)  point  of  z 

Load  Current 

pulse  width  (square  wave) 
open  circuit  voltage 
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